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Abstract. Fluid discs and tori around black holes are discussed within different 
approaches and with the emphasis on the role of disc gravity. First reviewed are 
the prospects of investigating the gravitational field of a black hole-disc system by 
analytical solutions of stationary, axially symmetric Einstein's equations. Then, 
more detailed considerations are focused to middle and outer parts of extended 
disc-like configurations where relativistic effects are small and the Newtonian 
description is adequate. 

Within general relativity, only a static case has been analysed in detail. 
Results are often very inspiring, however, simplifying assumptions must be 
imposed: ad hoc profiles of the disc density are commonly assumed and the effects 
of frame-dragging and completely lacking. Astrophysical discs (e.g. accretion 
discs in active galactic nuclei) typically extend far beyond the relativistic domain 
and are fairly diluted. However, self-gravity is still essential for their structure 
and evolution, as well as for their radiation emission and the impact on the 
environment around. For example, a nuclear star cluster in a galactic centre may 
bear various imprints of mutual star— disc interactions, which can be recognised in 
observational properties, such as the relation between the central mass and stellar 
velocity dispersion. 
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1. Introduction 

The interplay between gravitational attraction and non-gravitational forces produces 
systems with various shapes in the universe, including elongated structures, spheroids, 
discs, shells and rings. All of them participate in variety of processes and act on 
different scales, ranging from cosmological dimensions down to planetary sizes. Some 
configurations, spheroidal in particular, hold by themselves, while other are maintained 
by an imposed control of a central body or an external field. This paper deals with 
discs and tori — typical configurations of rapidly circulating fluids that are governed 
by a central object. Its gravitational attraction supports global stability of the 
configuration, at least below a certain critical density of the disc, while compactness 
of the centre suggests prospective observational methods. High density configurations 
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suffer from instabilities and their mass contributes to general-relativity effects jointly 
with the central mass. 

Nowadays it should be self-evident that black-hole discs are of indisputable 
relevance to different research communities, ranging from pure relativists to astro- 
physicists and observational astronomers. This diversification is reflected in varied 
approaches that people adopt to the problem. A mathematically minded theorist 
would take more interest in non-trivial (but still manageable) analytical solutions 
to Einstein's equations endowed with axial symmetry. An observer will want to 
actually see gas flows swirling in an accretion disc near a black hole (bh) or a compact 
star. Astronomers often invoke accretion discs to explain spectral features observed 
in galactic nuclei and in close (interacting) binary stars, although detailed physics 
of accretion flows still remains uncertain. Therefore, it is tempting to search for 
common intersections in such a broad subject. However, we will argue that serious 
distinctions cannot be neglected which prevent one to develop an all-encompassing 
scheme that would describe inherently different configurations. In particular, effects of 
general relativity (gr) have a chance to govern only the innermost regions of accretions 
discs and their role is quickly diminished with distance from the centre. Local (non- 
gravitational) physics plays crucial role and it often dominates over gravitation and 
governs observable properties of real systems (perhaps the most profound example 
arises if one attempts to compare accretion discs around stars with those around 
super-massive BHs). In order to see where limitations of any unified scheme stem 
from, it is useful to discuss both simplified and general analytical models and some 
particular astronomical objects in more detail. 

In this paper, two-dimensional models are discussed. Needles to say, this is 
already enormous restriction. Obtained solutions offer a convenient tool to elucidate 
physical mechanisms connected with non-spherical sources of strong gravitational 
field. Compact toroids represent a non-trivial model of such objects that require 
GR, although Newtonian approach provides a sufficiently accurate description in many 
cases. Interaction with the central field is especially intriguing if the latter is generated 
by an ultra-compact body. 

The paper starts with a summary of astrophysical motivations (section Our 
attention is focused on discs around super-massive black holes, because in this case the 
importance of self-gravity has been traditionally acknowledged, but we also mention 
compact binary stars, in which case it has been recognised relatively recently that 
self-gravity may also play a role, especially in the context of neutron star collisions, 
tidal disruption and subsequent emission of high-energy gamma-rays and, possibly, 
gravitational waves. In section |21 the problem of massive discs around black holes 
is discussed within GR. Whereas GR is mainly relevant for the inner parts of self- 
gravitating discs (within one hundred Schwarzschild radii, < lO^rs), the field of actual 
discs is equally important for their structure at larger distances from the centre where 
the Newtonian regime is adequate (section 0J|. Evidence for self-gravity in accretion 
discs is still rather tentative and incomplete; some possibilities and observational 
implications are explored in section [S] There we also indicate the potential role 
of other ingredients, such as a compact nuclear star cluster in which the accretion 
flow may be embedded. However, the reader should be aware that important parts 
of the whole debate remain untouched in the present paper. In particular, we do 
not discuss magneto-hydrodynamical processes, neither do we address non-axially 
symmetric modes of the disc flow. One should bear in mind that these processes are 
probably crucial for understanding the mechanism of viscosity as well as launching of 
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jets in astronomical objects. 

Disc structures are ubiquitous in such objects where the fluid circulates around 
and inflows onto a central compact body. The central mass, M,, can be considered 
as one of the model parameters and varied by many orders.^ Its value provides 
primary classification of cosmic objects. Physical characteristics of accretion discs 
scale roughly with M,. Indeed, accretion discs around super- massive black holes 
in active galactic nuclei and quasars share some properties with circumstellar discs 
in close binary systems (e.g., cataclysmic variable stars, X-ray binaries and micro- 
quasars). For example, both the discs around stellar- mass black holes and accretion 
discs around super-massive black holes may develop local and global instabilities which 
involve gravity. However, there are important distinctions between these two kinds of 
objects which prohibit any simplistic scaling. It turns out that galactic nuclear discs 
tend to be cooler and less dense compared to circumstellar discs, which translates to 
different opacity of the material, its ionization structure and, hence, different spectral 
properties. There has been also a very lively and still open discussion, whether the 
presence of a black hole in stellar-mass bh candidates (or, rather, the absence of any 
hard surface) can be revealed in the radiation properties of their discs. 

Environmental effects bring more distinctions between different objects with 
accretion discs. Unfortunately, at this point a bias is always inevitably introduced to 
our models by assumptions about turbulent transport of matter and momentum. Most 
solutions are model dependent. Neglecting this dependence can hardly be avoided 
when showing results of computations, nevertheless, we attempt to accentuate generic 
features when possible. 

We refer to monographs by Frank, King & Raine [113] and Kato, Fukue & 
Mineshige [170] for exposition of accretion theory, to Krolik [188] and Peterson [273] 
for discussion of accretion flows with the emphasis on active galactic nuclei, and to 
various review papers [7,73,174,209,266] summarizing recent advances as well as the 
problems encountered by traditional scenarios. 

2. Astrophysical motivations and some early works 

Accretion plays a crucial role in the processes of energy liberation and mass 
accumulation that take place in the core of astronomical objects. Disc-type accretion 
represents an important mode which is realized under suitable circumstances, defined 
by the global geometrical arrangements, local microphysics of the fluid medium and 
the initial conditions as well. Rotation is crucial and it puts discs, toroids and 
rings in the same category of centrifugally supported systems. By accretion discs 
one understands topologically toroidal fluid conflgurations orbiting around a central 
body and gradually drifting towards it, provided a mechanism to redistribute the 
angular momentum. The centre governs or at least substantially contributes to the 
gravitational field in which the medium evolves. The formation of the central body is 
presumed as being completed, although the accretion process gradually transfers mass 
and angular momentum onto the central body and thus contributes to its subsequent 
evolution. 

It is worth noting that the genesis of our own solar system with a star in the 
centre was imagined, by Pierre-Simon Laplace and Immanuel Kant, as emerging out 

^ We introduce Af, for the black-hole mass in order to maintain clear distinction from the disc 
mass. Both contributions add to the total mass of the system. Similar situation arises with angular 
momentum, as discussed later in the text. 
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of a turbulent, slowly spinning nebula with flattened, disc-like shape. ^ This was in 

the final years of the 18th century. The potential associated with a gravitating ring 
was examined in detail by F W Dyson [94]. Later, considerations associated with 
star and planet formation were the main driving force behind the initial development 
of accretion disc theory. Modern evidence for such configurations was reported as 
early as in 1930s by Baxandall [30] in the context of binary stars. In the attempt to 
interpret spectrograms of the famous /3 Lyrae binary system and in order to explain the 
origin of secondary lines, visibly displaced in their wavelength, "circulating currents 
in the outer envelope" were tentatively invoked together with a gaseous envelope that 
surrounds the larger component. Since then, interacting binary stars have been the 
prototypical systems for studying accretion discs. 

The treatment of self-gravitating discs was pursued in Ostriker's Newtonian 
equilibria of uniformly rotating, polytropic, slender rings [258]. Later, self-gravitating 
configurations with realistic equation of state and opacity were constructed and the 
basic formalism for self-gravitating bh discs was given in GR [27,46, 109]. In stellar 
binary systems, the gas in the disc is provided by a donor star via Roche-lobe 
overflow or stellar wind [180]. Accretion process helps to feed the central body. If 
conditions are suitable, the gas radiates and provides a way to observe and study the 
system. Historically, limiting cases were distinguished for methodological reasons and 
for the sake of simplification. When speaking of hydrodynamical non-sclf-gravitating 
accretion, the main assumptions concern stationarity, spatial symmetries that are 
imposed on the flow and on the gravitational field (spherical flows with negligible 
angular momentum, axially symmetric accretion onto a moving body, and disc-like 
planar accretion), and the sound speed profile (subsonic versus supersonic fiows). See, 
e.g. [42, 112, 146, 147] and references cited therein. It turns out that behaviour of truly 
three-dimensional unsteady flows is rather different and more complex [91,157]. 

Nearly spherical accretion occurs if the angular momentum content of accreted 
material is negligible (this case appears appropriate mainly for accretion onto an 
isolated body). In another limit, one ignores radial transport of the material and 
considers an axially symmetric toroidal configuration that preserves perfect rotation 
about the symmetry axis. Both these cases were treated in many works, assuming the 
test-fluid approximation (flow lines are determined by the central gravity and pressure 
gradients only). Toroidal topology of the fluid system introduces the main difference 
when compared with the physics of rotating stars (for the latter, see a text-book 
exposition of the subject by Tassoul [335]). On large scales, the accretion mechanism 
was proposed as the energy source in bright galactic nuclei in late 1960s (Lynden- 
Bell [214]). Since then the idea has been continuously developed and brought to a 
text-book level. 

Self-gravity has global consequences on the disc shape, namely, it inffuences the 
location of its inner and outer edges, as well as the disc geometrical thickness. The 
effect can be roughly parameterised by the ratio of the disc mass to the central 
BH mass, q = M^/M, [3, 185]. Wilson [363] speculated about the importance 
of self-gravity (especially in cataclysmic binary systems), which would have direct 
observational consequences, but this possibility was later rejected on the basis of 

^ "... We see a region of space extending from the centre of the Sun to unknown distances, contained 
between two planes not far from each other ..." (I Kant, Allgemeine Naturgeschichte und Theorie 
des Himmels, 1755). Independently, P S Laplace advanced and published a theory of Saturn's rings 
and the nebular hypothesis of the planetary system {Memoire sur la theorie de I'anneau de Satume, 
Mem. Acad. Sci., 1789; Exposition du systeme du monde, 1796). 
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more detailed arguments; it turns out that non-standard (massive) discs are ruled 
out in cataclysmic binaries [148]. Self-gravity is normally unimportant in binary star 
systems also because their discs have limited space to occupy. Under typical conditions 
(i.e. among well-known cosmic systems), the disc mass plays a role in the process 
of gradual fragmentation of proto-planctary nebulae, formation of planetary bodies, 
and their subsequent migration. (Recent observations, in particular with the Hubble 
Space Telescope, have demonstrated that discs are a by-product of star formation 
and also present during the birth of planets.) Also planetary rings are subject of 
self-gravitational instabilities — but we cannot tackle all these lively topics here 
(see [130,193,294,359]). There has been a recent interest in high-density (neutron) 
tori, in which case the disc own gravity must be important. This kind of massive discs 
could be formed in the course of binary star coalescence or as a product of collapse to a 
BH [135,243]. The question of stability and time evolution of these heavy, transient tori 
is essential. The idea was advanced in the framework of scenarios of failed supcrnovae 
and strongly magnetised tori [50,218,262,350,351], and an interesting possibility of 
intermittent accretion was hypothesised in the context of massive neutron tori that 
may be formed from a tidally disrupted compact star [371]. 

As mentioned above, gravitational field of heavy discs influences the global 
structure of the system, which in turn determines possible figures of equilibrium. 
Coming back to accretion discs in galactic nuclei, vertical structure of these discs is 
affected by their self-gravity especially in the middle and the outer regions (beyond 
a few hundreds rg, say) where the vertical component of the disc gravitational field 
exceeds the corresponding z component of the central field [316 318,329]. The main 
aspect of discs' non-negligible gravitational efi'ect on their own structure concerns the 
stability with respect to perturbations. Self-gravity was first examined in the context 
of accretion discs around super-massive black holes by Paczynski [260] , who proposed 
that local (Jeans-type) self-gravitational instabilities are invoked if the density exceeds 
a threshold value. Turbulent fragmentation may contribute to viscosity and to heating 
of the disc medium. Even the formation of stars may be induced in this way [186,319]. 
Collin & Zahn [76, 77] explored the possibility and the impact of star-formation in 
galactic-centre discs. Recently, also Sirko & Goodman [322] discussed consequences 
which the formation of massive stars in the disc may have on its heating and 
spectrum. Non-axisymmetric instabilities of accretion discs and tori have been widely 
investigated (see e.g. [39,129,267] and references cited therein). Their importance for 
modelling realistic accretion discs has been soon recognised, but we do not examine 
these types of instability here. If one considers strictly axisymmetric and stationary 
rotational motion, one finds idealized configurations, which provide useful insight into 
the problem of accretion discs; for example, the question of the role of mean stationary 
solutions. However, there is no doubt that these solutions must be rather different 
from real accretion flows; investigations of the latter has culminated in recent three- 
dimensional numerical simulations [91]. It is noteworthy here that we are allowed to 
treat both accretions discs in binary stars (with an accreting stellar-mass black hole) 
and in galactic nuclei, because common physical mechanisms operate in both cases. 
However, there are important distinctions which become particularly essential when 
realistic configurations are discussed. We will point out limitations of simple scaling 
between the two limiting cases later in the text, but great part of this paper concerns 
very simplified situations in which rough scaling according to the system size and mass 
is possible. 

In the next section we summarize techniques that have been employed to derive 
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Figure 1. Left; Contours of the Roche potential are plotted for two point 
masses ("stars") in their orbital plane. Here, the ratio of masses is assumed 
to be 0.1 (coordinates are scaled by separation between the two components; no 
disc is present). Right: The same type of graph as on the left side, but with a 
contribution of a thin planar disc taken into account as an additional source of the 
gravitational field. The plot shows the potential structure in the critical region 
near to the cusp, where the mass exchange takes place. See [363] for an original 
discussion of the binary Roche potential with a massive disc taken into account. 

spacetimes of gravitating discs around black holes within GR framework and then we 
proceed to possible astrophysical applications of massive discs. We start the discussion 
by briefly recalling a uncomplicated case of a light non-self-gravitating Newtonian disc 
that can be formed at a certain evolutionary stage of a binary star when the donor 
component dumps material on the other one. A text-book analysis of matter transfer 
employs the notion of Roche potential in which the actual gravitational field of a 
binary is approximated by the field of two massive points orbiting around a common 
centre of mass [45] (see the left panel of figure^. An accretion disc forms when 
the donor fills its Roche lobe and starts overflow near the Lagrange Li point. The 
gas transferred in this way possesses high angular momentum and creates a disc-like 
structure before being captured by the accreting body. The Roche model gives a 
basic classification of binary stars into three main categories according to the sizes of 
individual components: detached, semi-detached and contact systems. As mentioned 
above, their discs can usually be treated as non-gravitating. 

In contrast to this transparent scheme, a sufficiently massive disc could change the 
potential field near the cusp. In the right panel of figure ^ distortion of equipotential 
surfaces is depicted with the disc contribution to the gravitational field taken into 
account (here, the disc mass is comparable to that of the secondary star). It is relevant 
to remark that equilibria of relativistic tori also show the structure of critical lobes, 
which is similar to the Roche geometry, even if the fluid is non-self-gravitating and 
the only source of gravity is the central black hole (see section 13. 3|) . This structure 
facilitates accretion onto the central body. 

Equilibrium figures of critical (lobe-filling) tori represent a starting point for the 
discussion of accreting black holes in compact binary systems. Naturally, the Roche 
lobe geometry is applicable for the mass overflow onto a black hole as much as it 
is in other systems of close binaries, in which the mass is exchanged between the 
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components. However, topologically similar configurations of critical tori can arise 

also when accretion proceeds onto an isolated black hole. In this case the torus can 
develop a cusp and the overflow for a purely relativistic reason. The latter situation 
is particularly important for accretion onto super-massive BHs in galactic nuclei. 

An issue of global stability of such configurations must be raised: if a small amount 
of material is transferred through the inner cusp onto the central body, then a new, 
modified distribution of mass and of angular momentum is established. Geometrical 
shape of the critical equilibrium is changed accordingly. A question arises whether 
the new configuration is stable or whether, instead, it undergoes further accretion in a 
catastrophically increasing (runaway) rate. Abramowicz et al [2] claimed the presence 
of this kind of instability in pseudo-Newtonian massive tori, while Wilson [364] 
demonstrated that mass-less tori (i.e. test fiuid) around a Kerr bh are not run-away 
unstable. Answer to the problem of this instability depends on the radial distribution 
of angular momentum throughout the initial configuration (before the mass exchange 
takes place) and, in the relativistic case, also on bh angular momentum parameter. It 
was demonstrated that increasing bh rotation and a positive slope of the angular 
momentum distribution both have a stabilizing effect [6,86,110,213]. Increasing 
Md, on the other hand, acts against stability [228,229,253,254]. Comparisons [111] 
are not easy because different approaches have been employed in papers (numerical 
techniques vs. analytical calculations; proper GR treatments vs. approximations and 
pseudo-Newtonian models). The black- hole-disc systems are difficult to treat exactly, 
even in the test-fluid approximation. Severe simplifications have thus been imposed 
on theoretical models, ranging from presumed symmetries and boundary conditions to 
various assumptions about microphysics of the medium. Simplifying assumptions are 
the more unavoidable when the disc gravity (playing a role of the external source to the 
gravitational field) can not be neglected and when one intends to proceed analytically. 
Further restrictions are dictated by applicability of the results to real cosmic objects. 

3. Stationary axisymmetric fields in general relativity 

There are at least three aspects in which the system of a rotating black hole with 
an external source attracts relativists. First, due to the non-linearity of Einstein's 
equations, the field of a multi-body system is a traditional challenge where, in most 
cases, one does not manage with a simple superposition. On the first post-Newtonian 
level, the "celestial mechanics" of gravitationally interacting bodies can be kept linear 
(see Damour et al [87,88]), but in the strong-field regime the interaction may bring 
surprising features that have only been described in very few cases yet (for a two-body 
problem, the current "state of art" is the third post-Newtonian approximation [89]). 

The second point concerns manifestation of frame dragging due to the rotation of 
the sources. Contrary to the Newtonian treatment (that does not discriminate directly 
between static and stationary situations) , the relativistic field is determined not only 
by mass-energy configuration, but also by its motion within the bodies: the inertial 
space can loosely be imagined as a viscous fluid mixed by the sources. In today's 
"gravito-electromagnctic" language: gravity has not only an electric component, 
generated by the mass, but also a magnetic one, generated by mass currents [340]. 
Let us only refer to somewhat casual selection [68,97,163,176,216,227,271] for further 
details at this point. 

Needless to say, the third point is the presence of a BH itself. Near such an 
extreme body, the deviations from Newtonian theory become dominant; in particular. 
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the above mentioned impHcations of non-Hnearity and rotation reveal themselves 
prominently. 

The class of stationary axisymmetric solutions of Einstein equations is the 
appropriate framework for the attempts to include the gravitational effect of an 
"external" source in an exact analytical manner. At the same time, such spacetimes 
are of obvious astrophysical importance, as they describe the exterior of a body 
in equilibrium [158] (see [27,51,54,59,65,97,206,360] for thorough expositions on 
stationary and axially symmetric spacetimes in general relativity). 

In the Weyl-Lewis-Papapetrou coordinates {t, R, (j), z) of the cylindrical type, the 
stationary axisymmetric metric can be written as^ 

ds2 = -e^'^At' + B?B^e-^%dcl> - wdtf + e''^-^''{dE? + d^^), (l) 

where the unknown functions v, B, iv and A only depend on R and z; w is interpreted 

as the angular velocity of incrtial-frame dragging with respect to observers at rest 
at spatial infinity. The metric coefficients gab have invariant meaning, they can be 
expressed in terms of Killing fields k^^ = dx^ / dt and m'^ = dx^ / d(j): 

gtt = -e^'' + uj'^g<t><t> = hk\ (2) 

9t<t> = -^9<f<f = ^t"^'' (3) 

= R^B^e-^-" = m,m\ (4) 
The corresponding sub-determinant is det(gfa5) = —R^B"^, so 

g*t = _e-2-, 5*^ = -e-2v = -e-^'-w^ + (5) 

The Einstein equations read 

V-{RVB)=87rRB{TRR + T,,), (6) 

V • {R'^B^e-^''V(j) = -IGnBe'^^-^" r|, (7) 

V • {BViy) = iR^B^e-"^" {Vujf + 4nBe^^ (2T« + e'^"^ T^) , (8) 

where V and V- stand for the gradient and divergence in a flat three-dimensional space 

with coordinates {R,<j),z); thus VX = {X^r,0,X^^) and V • X = R-^[{RX"\r + 
{RX^-),z] in the axially symmetric case. Once B, oj and v are known, A can be 
integrated from equations 



(9) 



(10) 



BX^z — + R {b^rX^z + B^z^M. — B^Rz — 2Bv^rv^^ 
+ ^R^B^e-^'''uj^RUj,z = SttRBTr^, 

2BX^R — 2B^R + R 2B^rX^r — 2B^zX^z — B,rr + B^^z — 2B (iP^j^ — p^^) 

+\R^B^e-^- {lo% - ujI) = 87:RB {Trr - T,,) . 

The unknown metric functions arc subject to boundary conditions on the horizon (if 
there is one), on the symmetry axis and at spatial infinity. The metric must be regular 
on the horizon and on the axis, and in the case of an isolated source one also requires 
asymptotic flatness (the conditions are discussed in [27,54], for example). 

^ Geometrized units are used in which c = G = 1; the signature of the metric tensor p^i/ is ( — !-++)• 
Greek indices run from to 3 and Latin indices (i, j, ■ ■ ■) run from 1 to 3; indices from the beginning 
of the Latin alphabet (a, 6, ... ) will represent cyclic coordinates t and 4>. Partial differentiation is 
denoted by 9 or by a subscript comma, covariant derivative is denoted by V. 
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Outside of the sources, where T^^ — 0, the simplest solution of equation (j^ with a 
satisfactory asymptotic behaviour (equation H2t)|) below) is i? = 1. The two remaining 
field equations lEJ)~l|Hl) reduce to 

V • (R'^e-^'^Vu;) = 0, (11) 

VV= ii?2e-4-(Vcc;)2 (12) 

and the relations iPl- lfTUI) yield 

A,;,, = R ~ vl) - ii?3e-4^(a;2^ - u^^,), (13) 

\,^2Rv^RV^,-\R''e-^''uo^RUo^,. (14) 
Equations lfTT)l - (fT^ are often presented in the form of the Ernst equation [98] 

- 5tt W = {VSf (15) 
for a complex (Ernst) potential S = —gtt + iV" whose imaginary part ip is given by 

R'ip.z = gttgt4>,R- gtt.,Rgt4„ R'4',r^ gtt,z9t<i,- 9tt9t4>.z- (16) 

A is found by a line integration as above. 

3.1. Essential features of the stationary axisymmetric spacetimes 

A black-hole horizon is a null hypersurface below which the spacetime is dynamical 
(see e.g. [54,55,117] for thorough accounts). In Weyl-Lewis-Papapetrou coordinates, 
the horizon is located where det{gab) = —R?B^ vanishes. In the vacuum case we can 
put i3 = 1, so the horizon lies on the axis then, i? = 0. The BH interior thus has to be 
studied in different coordinates, e.g. in spheroidal coordinates of the Boyer-Lindquist 
type {t, r, 0, 0), introduced by the transformation 

i? = VAsin6l, z^{r~M,)y, (17) 

where A = (r — Af,)^ — fc^ and y = cos6', M, being a scale parameter (it represents 
the black-hole mass) and (±)fc determining where the horizon reaches up along the z 
axis. The horizon is now given by A = 0, i.e. r = M, + k = r^ (hence, z = fcy; there 
exist more black- hole horizons in general, but we mean the outermost one here). 

Important parameters of the horizon are its surface area A, surface gravity Kh 
and the angular velocity ujn relative to infinity. 



\l {gee g4,4>)ii^0, Kh = [Vc'^Ih, lu^ ^ uj^r ^ r^i); (18) 







Kh and iOii are constant all over the horizon. 

In order to learn a "true shape" of the bh, not distorted by coordinates, 
the horizon must be represented properly as a two-dimensional surface in a three- 
dimensional Euclidean space. One can e.g. follow [323]. The two-dimensional metric 
is first rewritten as 



H 2 A 



h-'iy)dy' + h{y)d(b' , (19) 



where h{y) = A7rA~^{g^cf,)-ii. An isometric embedding of the two-surface (y, 0) in E'^ 
with coordinates {X, Y, Z) is then given by 

y I 

X^^Vhcos^, Y=^Vhsinc^, Z=^JJl[l-h^yy. (20) 

* 
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2 2^ . rr. . (24) 



If the horizon's Gaussian curvature 

"^^^ -8n^A-^h^yy (21) 

is negative in some place, the global embedding is impossible. 

On the symmetry axis {R — 0), the regularity condition requires e^ = B, and 
therefore the metric adopts the form ds^ = —e^^dt^ + B^e^^'^dz^ here. Simplification 
of the formulae can also be expected in the equatorial plane {z = 0) if the spacetime is 
reflectionally symmetric with respect to it. This is usually the case for astrophysically 
motivated considerations. 

In the static case uj — 0. If + — (which is fulfilled with zero pressure), 
the metric acquires the Weyl canonical form 

= -e^-dt^ + i?2e-2-d02 + e^^-^-" (di?^ + dz^) . (22) 

Equations l(B|)- (|10(l reduce themselves to Poisson's equation 

= + ly^RR + v^,_, ^ Aire^^-^'iT^ - T,*) (23) 

for ly, and to relations 

\.z -^Rv.RV^z = 8ttR Trz , 

X.R - R {u% ~ ul) = 47ri? {Trr ~ T, 

for A. The linearity of Laplace equation allows for simple superposition of vacuum 
Weyl fields. 

From the astrophysical point of view, a given analytical solution is problematic if 
it contains features which should not be present in a realistic source, such as (i) physical 
singularities on or above the horizon, or (ii) bad asymptotic behaviour. The existence 
of a singularity follows — according to its type [69, 95] — e.g. from the divergence 
of scalars constructed from the metric tensor or from its derivatives, or from the 
divergence of physical (tetrad) components of the Riemann tensor. The Kretschmann 
invariant R^upaR^'^'^'^ is usually checked first. However, the spacetimes are known 
which do contain singularities, although all of their curvature invariants vanish. It is 
thus difficult to verify the non-existence of singularities; the only proof of regularity 
at a given point is to find a coordinate system in which the metric is smooth enough 
locally. 

Astrophysically relevant is the region of outer communications outside the 
horizon, but a relativist is interested in the black hole interior as well, because even 
there it is possible to perform physical observations. In the Weyl-Lewis-Papapetrou 
coordinates the metric is not defined below the horizon (this would correspond to 
imaginary radius R) and must be extended there. One asks about how is the bh 
interior (and the singularity in particular) influenced by the presence of external 
sources. The answer is unknown even for quite simple superpositions. 

Far away from the horizon, one inquires for the asymptotic behaviour of the field. 
That of an isolated stationary source falls to zero in a specific manner — the spacetime 
is said to be asymptotically flat (e.g. chapter 11 in R Wald's book [358], or section 3 
in J Ehler's Festschrift [32]). In case of the metric (Ql, it must hold, for r oo, 

total mass ^, _2\ 

)> (25) 



B = 1 + &{r-^), (26) 

2 • (total angular momentum) , , ^ 

UJ - — ^ + e{r-^), (27) 
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A = ^(r-2). (28) 

Although a suitable combination of solenoidal motions can satisfy the assumption 
of stationarity and axial symmetry, it is natural to assume that the elements of the 
source follow the simplest type of motion — that along spatially circular orbits {R = 
const, z = const) with steady azimuthal angular velocity Q = d(j)/dt. Such a motion 
follows the symmetries, so an observer on a circular orbit experiences time- independent 
field around. In order that the corresponding four- velocity = -I- ilm'^) points 

inside the light cone, (m')^^ = —gtt — '^^9t<p — ^'^g<p<p = e^''(l — v"^) must be positive 
(we denoted v = RBe~^'^{fl — u) the linear speed with respect to the local zero- 
angular-momentum observer). The latter condition is fulfilled if Q falls within the 
range 

e" 

Omin = UJ - — — <D, <UJ + — — = rimax- (29) 

On the horizon (where e^ = 0), the permitted range narrows down to a unique value 
u)h- In the above-given equations, the four- velocity covariant components 

7 = -ut = -u*{gtt + ^gt4>) = u*e^'' + ^i-, 
I =u^ = u*500(f2 - w) 

stand for the specific energy and the specific azimuthal angular momentum with 
respect to an observer at spatial infinity. 

In case of the circular orbit, the four-acceleration = Uf_i-„u'^ can be written as 

a^--25a/3,M" 2e^" (1 - v^) ' ^ ^ 

It has at most two non-zero components, a/j and o,^; in general, is zero in the 
equatorial plane and non-zero elsewhere, while aji vanishes if 



Q = = -9t,p,R ± \/ {gt4>.Ry - gtt,R g<t><t>,R ^^2) 
g<t><t>,R 

where the upper/lower sign corresponds to a co-rotating/counter-rotating orbit (these 

terms may be somewhat misleading in spacctimes with large angular momentum). 
Among equatorial geodesies, three particular cases are important: the photon, the 
marginally bound and the marginally stable circular orbits. The photon orbits 
delimit the regions where the circular motion is time-like; they arc given by equalities 
J7± = f2max. The marginally bound orbits demarcate the regions where particles on 
circular orbits have lower energy than necessary for the existence at spatial infinity; 
the limiting case is given by ^{Q, = Q.±) = 1. The marginally stable orbits also bound 
the regions where circular motion is stable with respect to small perturbations acting 
within the orbital plane. Rayleigh's criterion of linear stability gives the location of 
the marginally stable orbits where = ^±)]^r = 0. These orbits play a key role 
in the theory of accretion discs: the matter is swept away from unstable sectors. In 
particular, the innermost marginally stable orbit should represent the disc inner rim. 
However, perturbations perpendicular to the orbital plane can also be important.** For 
the above-mentioned exact solutions of Einstein's equations, the frequencies of vertical 
oscillations were given explicitly [311]. 

■* It has been discussed quite recently that (non-self-gravitating) accretion discs are prone to warping 
due to extraneous irradiation [281]. This effect is, naturally, modified when self-gravitation is taken 
into account [269]. 



Topical review 



12 



A peculiar feature of rotating fields are the dragging effects. In accordance with 
Mach's ideas of relativity of motion and inertia (and his interpretation of Newton's 
bucket experiment in particular), Lense & Thirring showed, in the early times of 
general relativity, that dragging is present in Einstein's theory. Machian inspiration 
remains alive within contemporary relativity — see, for example, the last-decade 
references [68,163,176,215,216] or the Tiibingen conference [25] and references therein. 
The effect has obvious analogy in electromagnetism, where (electric) currents generate 
the magnetic component of the field. An extreme implication of dragging is the 
occurrence of the ergosphere in the vicinity of ultra-compact rotating objects. In 
this region f2niin > 0, thus the stationary observer cannot remain static (i.e. at rest 
relative to infinity), albeit he still withstands the radial attraction. The static-limit 
surface, given by fimin = and hence by gtt = 0, also represents the set of points 
from where signals emitted by static observers (fl = 0) reach infinity with an infinite 
redshift. Several processes working in the ergosphere have been suggested, by means 
of which the rotational energy of the black hole could be extracted without violating 
the second law of bh dynamics [283]. 

With the metric Q, 5tt = corresponds to RBoj = e^^ which can only be 
satisfied by i? > 0, so the static limit really lies outside the horizon;'^ in particular, 
i? = implies e^" = 0, so the static limit touches the horizon at the axis. 



3.2. Sources 

In order to complete the relativistic solution, one must also describe the interior of 
sources — where T^^ ^ 0. It is difficult to find a realistic interior solution, especially if 
it has to match smoothly the vacuum exterior and if we restrict ourselves to analytical 
approaches. The problem can be substantially simplified by assuming that the source 
is infinitely thin: the solution is then vacuum-type everywhere, with the energy- 

— 1/2 

momentum tensor T^i, — g^z Sfj_i,5{z) found from the discontinuity of the normal 
field across the source as in electrodynamics (the appropriate covariant method is 
known as the Israel's formalism [29,159]). 

The case of a thin equatorial layer (a disc) in a stationary axisynimetric spacetime 
has notably been studied by [126, 197]. The disc lying on z = const has zero radial 
pressure inside {Sua = 0). At the same time, the vertical pressure vanishes {S^z = 0) 
if the disc mid-plane coincides with the equatorial plane (z = const = 0) of the 
reflectionally symmetric spacetime. In the Weyl-Lewis-Papapetrou coordinates, the 
non-zero S^i, components read [197] 

g _ ^/9RR f gab \ 

Using and ((nil , this yields 



(33) 



4//,,(l - Riy^R) - R'e-^-'u^.iio - Ruj,r) , (34) 



Stt 

Si = -^i?2e-4v., (35) 
^ Stt 

sf = [Alov,, - (1 + R^e-^-uj^)uj.,] , (36) 

87r 

^ In coordinates which also describe the BH interior, two horizons and two static limits are found in 
general. We mean the outer horizon and the outer static limit everywhere. 
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'iv,RV.z- Re-'^''uj.^{Lj + Roj.r) . (37) 



The total mass and angular momentum of the disc are fixed by the Komar integrals 
that, respectively, lead to 



Md g'^gtc.RdR = -i Z ^ ) di?, (38) 



OO ^2 



^ ^ 9''9c4,,zR di? = i y ^ c^,, dR; (39) 

b 

R = b > is the position of the disc inner rim. All the z-derivatives are understood 
to be calculated in the limit z ^ 0~^. 

If a BH is present in the disc's centre, its contribution to the total mass and 
angular momentum need to be included. For a given type of spacetime, the aggregate 
values of the parameters can be written as sums, Af, + Md and J, + Jd (cf. [55], 
chapter 6.6.1). The black-hole parts can be calculated as Komar integrals, or they can 
be inferred by subtracting Md, Jd from the totals, derived from the metric asymptotes 
H25I27|I . M, and J, are related through Smarr's formula,^ 

M, = 2whJ. + — = 2whJ, + fc . (40) 
Further, if the term 

D ^ [st - Siy + ASlSf = {4^ J - R'e-^-'u;^^) (41) 

is non-negative, then a stationary observer exists with respect to whom the energy- 
momentum tensor assumes a diagonal ( "isotropic" ) form 

S^" = wu^^^ul,, + Pv^^Xso, (42) 

where = ^'^^(l, 0, fiiso, 0) and vt"^^ = -R~^{£iso,0,Jiso,0) are, respectively, 

observer's four-velocity and the unit basis vector in the (/)-direction. Angular velocity 
of this observer is 

".-^(«:--;-v^^"-IS:-\/(l&rS^ 

7iso and ^iso denote the corresponding specific energy and angular momentum at 
infinity.'' The observer measures the surface density 

2 Ctt _ f2 q4>(t> 

"'iso V /ISO ^ ^ ^iso^iso ) 

and tangential pressure 

p_ R^uUs<f^-nl^s'*) 

[ISO ^ ^ ^ISO^ISO 

^ See Carter [55], equation (6.297), for derivation and discussion. For the original formulation in 
Kerr spacetime, see [324]. 

The case of D < represents discs with non-zero heat flow. Then the energy-momentum tensor 
must be written in a more general form, S'^'^ = wuf^ uY + Pv^ vV + K(uf^ vV + uV ), with 

' ISO ISO ISO ISO ^ ISO ISO ISO ISO'' 

K = s/-D/2 and Qi^a = (25^-1(5^ - 5*) = ui - 2R-'^u)-^e'^''v^^. See [126]. 
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li w > P > 0, the energy- momentum tensor can represent two equal streams of 
particles moving on (accelerated) circular orbits in opposite directions at the same 

speed ij P/w. This "speed of sound" was shown [126] to be equal to the geometric 
mean of the local prograde and retrograde circular geodesic speeds v± as measured by 
the observer u^^: 

\v+V-\^P/w, (46) 

where 

(Wiso)t + (Wiso)0f^± . 
(Uiso)t + (Uiso)0"± 

r2± are the angular velocities of the prograde and retrograde circular geodesies 
Although counter-rotating relativistic streams are hardly applicable to any realistic 
astrophysical situation, they can be very illuminating especially in the discussion of 
physical sources of the Kerr metric and other spacetimes.* 

Further physical requirements can be raised, namely, the weak energy condition 
{w > 0, P > —w), the dominant energy condition {w > 0, |P| < w) and the non- 
negativity of pressure (P > 0). Their combination gives w > P > 0. The requirements 
are very simple in the static case: if I'^ziz = 0+) > 0, then P > is equivalent to 
'^,R > which is satisfied below the Lagrangian point of zero field (where ly.n — 
and Vj^ =0). In regions with tension (P < 0) the particles would be attracted more 
by the outer parts of the disc than by the central bh, thus no circular geodesies exist 
(then hoop stresses have to be employed to interpret the disc matter); this typically 
happens in discs with much matter on large radii. The other condition, w > P, is 
equivalent to < 1. To summarize, the above constraints are in fact ensured by the 
obvious condition < v± < 1 which can be written explicitly as 1 > 1 — Riy^B. > 1/2- 



3.3. Solving Einstein's equations for a black hole with matter around 

Relativistic spacetimes are being found in three ways: (i) by numerical solution of 
Einstein's equations, (ii) by perturbation of previously known spacetimes, and (iii) by 
exact analytical solution of Einstein's equations. Let us mention several results of the 
above approaches on stationary axisymmetric spacetimes that describe rotating black 
holes with additional matter. 

Teukolsky [338] lists today's computational resources and the appropriate, 
hyperbolic formulation of the field equations as the main "reasons why we are on 
the verge of important advances in the computer solution of Einstein's equations" . 
Whereas analytical prospects are restricted, the present-day computers can handle 
almost any situation. However, covering a sufficiently representative time-span is 
still a problem. Furthermore, it is often impossible to recognise whether a given 
numerical solution represents a typical situation or just a marginal case. Due to this 

® Bicak &; Ledvinka [35] have constructed pressure-less counter-rotating discs in the Kerr geometry, 
employing Israel's technique [29,159]. Analogous approach was used by other authors in a different 
context of exact solutions; see e.g. Abramowicz, Arkuszewski & Muchotrzeb [1], Turakulov [347], and 
Pichon & Lynden-Bell [275], or a particularly lucid paper by Bicak, Lynden-Bell & Katz [37]. These 
works may appear rather academic for people who are interested mainly in astrophysical applications 
of gravitating discs. However, a very pertinent question is addressed in these papers, namely, under 
what conditions a material disc could represent a source for the Kerr metric. Very recently, the 
"displace, cut and reflect" method has been extended to generate static thin discs with halos [355] 
and static thick discs [127]. 
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lack of generality, it is difficult to discern, analyse and interpret different classes of 
solutions across boundless ranges of possibilities. Nevertheless, numerical solutions 
provide explicit examples of spacetimes and processes that might otherwise remain 
only conjectural [198]. 

Numerical spacetimes with a rotating bh surrounded by a stationary 
axisymmetric source were already constructed by Lanza [191] (a hole with a thin finite 
equatorial disc) and by Nishida & Eriguchi [252] (a hole with a thick toroid); recently, 
a new ( "multidomain spectral") method has been utilized to spheroidal as well as 
toroidal fluid configurations [16], with the possibility of inclusion of a central object in 
future. To mention just one particular point, contrary to a common experience (also 
gained by [252]) that the horizons inflate towards the external sources, [191] ended 
with a prolate horizon (stretched along the rotation axis) in certain cases (when the 
disc was strongly counter-rotating with respect to the hole). It would be an interesting 
consequence of the interplay between dragging from the hole and from the disc if it 
were confirmed that this can indeed happen, even though under extreme circumstances 
only. 

A great deal of literature and several formalisms have been devoted to 
perturbations of BH spacetimes. Will [360] gave an explicit result for the perturbation 
of the Schwarzschild metric by a (rotating) axisymmetric weakly gravitating thin 
equatorial ring. This direct approach is however not practicable for a rotating 
hole and/or for an extended source with pressure (not mentioning more complicated 
situations). 

The rotating case (specifically, the algebraically special vacuum case) was 
discussed notably by Teukolsky [337] who succeeded in separating the decoupled 
equations for the first order perturbations of a Kerr BH into the second-order 
ordinary differential equations. By solving the "gravitational" Teukolsky equation, the 
perturbative deformation of the Kerr horizon was calculated by [90]. The stationary 
axisymmetric Green's function of the Teukolsky equation was provided by [210]. The 
most comprehensive expositions of the first-order bh perturbation theory were given 
by Chandrasekhar [60, 61] in the Newman-Penrose formalism. 

In a related approach to the first-order perturbations of rotating fields [66, 171, 
343,345], the perturbation components arc expressed in terms of a single (Debye) 
potential that obeys a wave- like equation (adjoint of the Teukolsky master equation). 
The problem of perturbation of the Kerr horizon was treated this way by Chrzanowski 
[67]. Yet another approach, not restricted to algebraically special backgrounds, was 
advanced in [301]. Note that all these perturbative methods deal with curvature 
components, as opposed to the earlier metric formalisms of Regge & Wheeler, Zerilli 
and Moncrief. 

In the charged case, the perturbation problem introduces coupling between 
gravitational and electromagnetic quantities. Conversion between gravitational and 
electromagnetic waves takes place. In the case of Reissner-Nordstrom BH, these 
interacting perturbations were studied by Bicak [33] and, more recently, by Torres 
del Castillo et al [344] . A gauge invariant derivation of the basic equations was given 
by Fernandes & Lun in Schwarzschild space-time [106], while [107] discussed a rotating 
analogy and provided its link with the Teukolsky equation. 

The second-order perturbation theory has also been under development since the 
1970s. For Schwarzschild black holes, it is surveyed in Gleiser et al [122], while the 
rotating case is tackled in Campanelli & Lousto [52]. 

It should be remarked that solutions describing stationary sources around black 
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holes are not the primary aim of the perturbation strategy. Historically, attention 
has rather been devoted to the stability of black-hole solutions, to the relaxation of 
black holes into a stationary no-hair state, to non-stationary processes of astrophysical 
significance (such as the behaviour of weakly gravitating particles or waves in BH 
backgrounds) and to gravitational waves (sec Thornc [339] for a survey). 

One can add that the approximate (though non-perturbative) method for 
superposing two (rotating) black holes was proposed by Krivan & Price [187] to fix 
initial data for both analytic and numerical computations of bh collisions. 

A perturbative approach is adequate in situations where the outside source has 
only a very small effect on the field of the main body. Whenever this is not the case, 
the superpositions have to be described by solutions of the hill, non-linear theory. At 
the end of Bleak's survey [34], one is encouraged "not to cease in embarking upon 
journeys for finding them, and perhaps even more importantly, for revealing new roles 
of solutions already known" , because "Is there another so explicit way of how to learn 
more about the rich possibilities embodied in Einstein's field equations?" More recent 
reviews of vacuum fields can be found in Bonnor et al [43,44]. For further details, 
sec also the canonical monograph on exact solutions to Einstein's field equations by 
Stephani et al [328]. 

Until now, no explicit exact solutions describing the system of a rotating centre 
with an additional axisymmetric ring, disc or torus have been found. Nevertheless, 
many stationary axisymmetric solutions of the (electro-)vacuum Einstein equations 
are known that do generalise those containing only isolated black holes. These were 
mostly obtained by indirect methods known as "generating techniques" . Two major 
approaches, developed by the end of the 1970's — the group-theoretic techniques 
and the soliton-theoretic (or inverse-scattering) techniques, work for spacetimes with 
two commuting symmetries. (Sec Kordas [181] for a review and [80-82, 144, 219] 
for more detailed analysis and interrelations between different formulations, e.g. 
those of Kinnersley & Chitre; Maison; Belinskii & Zakharov; Harrison; Hoenselaers, 
Kinnersley & Xanthopoulos; Hauscr & Ernst; Neugebaucr; Kramer & Neugcbauer; 
or Alekseev.) Other related methods have been proposed more recently, e.g. the 
simplification of the Hauser-Ernst integral equation by Sibgatullin [320], "monodromy 
data transform" by Alekseev [9] (also [179]), static gravitational multipolcs [131] 
(cf. [284]) and their superposition with stationary fields [132] by Gutsunaev & Manko, 
linear transformation by Quevedo [284] or "finite-gap" (algebraic-geometric) solutions 
by Korotkin & Matvecv [182]. Most of them have been compared in reference [10]. 
(For other approaches, see e.g. [71,334] or the results of Nakamura and Kyriakopoulos, 
referred to and worked out by [336]; cf. also [352].) 

The crucial point of the soliton generating techniques is a solution of two linear 
differential equations (Lax pair) the integrability conditions of which are exactly the 
Einstein equations (namely the Ernst equation). The linear problem can be tackled 
in order to generate new solutions from the known ones: given some ("seed") metric 
with two Killing vectors, it yields a different metric of the given type (the procedure 
is often called Backlund transformation in analogy with the technique used for the 
KdV equation). In such a manner, many known spacetimes were reproduced, but also 
broad families of new solutions were provided characterised by arbitrarily large sets 
of free parameters. Only a very restricted number of them have been given a clear 
physical interpretation, however. Though several of these results perhaps represent a 
rotating bh in an "external" gravitational field (e.g. [47,62,63,220,221,285]), none of 
the latter has yet been firmly linked with a concrete realistic body such as ring, disc 
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or torus (but cf. section IV in Letelier fc Oliveira [204]). Namely, when performing the 
Backlund transformation, it is difficult to require specific physical properties of the 
spacetime being generated — one rather tries what will come out. See [284] for a hint 
of how to overcome this. More recently, Letelier [203] interpreted a solution obtained 
from a Weyl seed by the inverse-scattering transformation with an even number N 
of solitons as a superposition of this Weyl solution — i.e. of a given set of static 
multipoles — with N/2 Kerr(-NUT) solutions, represented by N/2 rotating axial bars 
in Weyl- type coordinates. 

As an example, let us give the attempt [372] which was "cultivated" from a 
general Weyl metric by a real-two-soliton version of the Belinskii-Zakharov method. 
The resulting metric was written as a generalisation of the Kerr ( -NUT) solution and 
many of its properties were proven satisfactory. Choosing the inverted first Morgan- 
Morgan disc (see below) as a seed, various characteristics of the superposition were 
plotted [307] against the rotation parameter of the central hole, the mass of the disc 
and its inner radius. The characteristics simplify very much in the limit when the 
BH becomes extreme. One of the interesting features of this limit is the vanishing of 
the "external" gravitational field on the horizon [305] . This effect of expulsion of the 
external (stationary axisymmetric) fields from rotating (and/or charged) black holes, 
analogous to the Meissner effect in (super) conductors, was observed in magnetic fields 
before (see the recent survey [36]). However, a serious problem has been found in the 
equatorial plane: the black- hole horizon has a sharp narrowing there. This is caused 
by the appearance of a supporting surface between the hole and the source outside 
it [307]. Indeed, the possible problems with physical interpretation of the solution (e.g. 
negative density or pressure or matter moving at a superluminal speed) also involve 
the occurrence of supporting singular surfaces or lines ( "struts" ) that indicate that a 
given system of bodies cannot remain in stationary (or static) equilibrium. 

The struts may equally well "spoil" a static superposition, though this case almost 
reduces just to adding the two potentials, v = ui+ 1/2. The singularity may appear in 
calculating the second metric function, 

R,z 

\= j R (!/2j _ jy2J ^ 2v^R v^^ dz . (48) 

axis 

Consequently, even the static axisymmetric case has only been afforded a few realistic 

superpositions. Let us refer to two of them involving a Schwarzschild BH: the 
superposition with an infinite annular thin disc, obtained by Lemos & Letelier [201] 
by inversion of the first counter-rotating finite disc of Morgan & Morgan [241] and 
the one with an inverted isochrone thin disc, obtained by Klein [177]. The Morgan- 
Morgan family of static axisymmetric finite thin counter-rotating discs [241] is given 
by Newtonian surface density 



Since the zeroth member is clearly singular at the rim, most attention was devoted to 
the first (m = 1) member. It appeared to be a prototype of a simple and physically 
meaningful Weyl field and it was also considered as a limit case of a more general 
class of stationary discs [14,115]. The corresponding annular case, obtained by the 
inversion R jji^^'^ , z p^+z^ ' '^^^ shown to really allow for physically satisfactory 
situations [201,310,312,313,370]. There is nevertheless a problem at the rim of 
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these discs, already suspected by checking the radial derivatives of H49|l : the m-th 
member of the family has infinite m-th and higher derivatives of the density at the 
outer rim. Calculation of the Kretschmann curvature invariant reveals that the rim 
is really singular and that the singularity is inherited by the inverted discs [304,370]. 
The superpositions of a Schwarzschild bh with (inverted) "higher" Morgan-Morgan 
counter-rotating discs were studied in [308]. Comparison of the results obtained for 
the first ten of them showed that the properties of the disc are very sensitive to the 
local density profile. It will be desirable to also consider different classes of annular 
thin discs in superpositions, preferably of those not having any such drawback at the 
rim. 

We have already explored one different family of annular disc solutions, namely, 
the case when density is a power-law in Weyl radius, 



However, these discs also turned out not to be regular at the rim; (m -I- l)-st radial 
derivative of their potential diverges, although all the radial derivatives of density 
are finite there, and the first m — 1 derivatives even vanish. In a different context 
of galactic stellar discs, an interesting approach has been developed by SibgatuUin, 
Garcia & Manko [321] who examined a new method of reconstruction of the surface 
density distribution of a finite size disc with a central bh. 

Allowing for self-gravitation of the matter around the horizon, one has to reduce 
the generality of the problem elsewhere. Whereas test discs are treated within 
stationary axisymmetric settings, analytical solutions with exact inclusion of the 
ambient material have only been studied in more detail in static cases yet. This 
is of course a serious limitation, because the collapsed cores of galaxies as well as 
ultra-compact components in binaries are likely to rotate, perhaps very rapidly, which 
means that inertial dragging should be important there in addition to the effects of 
pure scalar attraction. Moreover, the composite spacetimes presented in the literature 
are constructed with discs whose field can be expressed in a tractable form rather 
than with those which really follow from some model of accretion. However, some 
superpositions — such as e.g. those involving thin annular discs obtained by inversion 
of the Morgan-Morgan solutions [201,310-313,370] — were shown to possess physically 
acceptable properties within an astrophysically plausible range of parameters. As such, 
they can indicate some of the effects that could occur as a consequence of non- negligible 
gravitational contribution of the disc (see [306] for a more thorough survey). 

Let us return to the linear problem (Lax pair of linear equations) associated 
with the Ernst equation. In order to control the physical properties of the solution 
being constructed, it is desirable to translate the physical boundary conditions of the 
Ernst equation into the language of quantities which appear in the linear problem and 
then solve the latter (rather than to apply a "random" Backlund transformation to 
some metric [202]). This leads to Riemann-Hilbert problems known from the theory 
of completely integrable differential equations [17,179,249]. Tackling the stationary 
axisymmetric boundary value problem with the help of the linear system, the field of 
two physically relevant types of sources has been discussed: that of black holes and 
that of finite thin discs ( [14,115,232,250,251] and references therein). The hope that 
the above methods could be used to describe superpositions of both (e.g. [250,251]) 
has recently been supported by Klein [178] who presented an exact space-time of a 




(50) 
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regular black hole surrounded by an infinite annular disc as a subclass of Korotkin's 
solutions to the Ernst equation, obtained by Riemann-surface techniques [183]. The 
resulting metric is given in terms of theta-functions, being asymptotically flat under 
suitable choice of parameters. Klein's solution seems to involve cases with physically 
acceptable source disc. 

The limit case should be briefly mentioned when gravitation of the external 
matter is ignored and the central rotating bh fully determines the field. Then the 
spacetime geometry is described by Kerr metric and the equilibrium configurations 
can be found analytically [26,109], at least if the fiuid obeys a simple equation of state 
(such as a barytropic relation) and remains in dynamical equilibrium with purely 
rotational motion (four- velocity u — u'^Bcj, + u^et). This case is over-simplified from 
the astrophysical point of view, but it captures several crucial relativistic effects which 
do not occur in the Newtonian limit — namely, the possibility of fluid overflow across 
the inner edge, which acts, effectively, as the Lagrange Li point in a binary system. 

The relativistic Euler equation can be integrated in terms of specific enthalpy, 
W = — /(e + P)^^AP — \n\ut\ — F{1), which can be further expressed as a 
function of /(i?, z), the profile of angular momentum density [4, 5, 12, 185]. The only 
unknown function that needs to be specified is the rotation law, either in terms 
of I — —u^/ut (which is a conserved quantity), or angular velocity ft = u'^/u*. 
Rotation can be differential but the choice is not completely arbitrary because, 
according to Rayleigh's criterion for linearly stable configurations, / must be increasing 
monotonically outwards from the rotation axis (a borderline case of constant I was 
discussed in detail). W plays the role of an effective potential for inertial forces and 
determines equilibrium configurations in which the equipotential and isobaric surfaces 
coincide, so this way constant pressure surfaces can be found. For example, in a special 
case of marginally stable configuration with / = const, one finds W{R, z) =h\\ut\. 

The fluid can establish a stationary configuration provided that the VF-isosurface 
is closed (see figure ISJ, otherwise overfiow and subsequent accretion take place. Let 
us recall that equilibrium configurations with their elongated cross-section (bounded 
by the critical, self-crossing lobe) resemble the shape of Roche lobes. The possibility 
of fluid transport across the cusp occurs when the critical lobe is filled; accretion can 
thus proceed without a necessary action of viscosity. However, the two situations are 
of very different nature. In the former case, it is the combined gravitational effect of 
the binary together with (Newtonian) centrifugal force which define the lobes, while, 
in the latter case, the only source of gravity is the central body with relativistic effects 
playing a crucial role. 

3.4^. Gravitational effects of black-hole discs 

Compared with the vacuum case (with a bh but without a disc) , modifications of the 
observable characteristics can be linked with the presence of self-gravitating matter. 
Various approaches have been adopted in order to construct solutions and discuss 
motion in the gravitational field of black-hole discs. One can start, e.g. , by constructing 
equilibrium structures of self-gravitating polytropic thick discs [40]. Formalism was 
developed also for perturbative solution of a (slowly rotating, weakly gravitating, 
geometrically thin) ring around a Schwarzschild BH [360,361]. Furthermore, equatorial 
motion was explored in pseudo-Newtonian studies [172] of a rotating bh surrounded 
by a massive thin equatorial ring, and in GR analysis of the corresponding static 
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Figure 2. The equipotential surfaces of effective potential W{x, z) are shown in a 
meridional section {x, z) of the Schwarzschild metric. Left: Contours ofW = const 
are plotted with solid lines. The critical self-crossing contour, W = —0.065, is 

indicated with a dashed line. Equilibrium is possible if the fluid (in pure rotational 
motion) fills one of the closed surfaces. Notice that these configurations can extend 
far from equatorial plane. The black hole horizon is indicated by a circle of radius 
7^ = 2M around the origin. Right: A detail is shown of the structure of contours 
near the cusp. Fluid overflow provides the possibility of relativistic accretion 
across the cusp without the need for viscosity. This effect also stands behind the 
origin of the runaway instability [5, 185]. 



case [57].^ The accretion-disc self- induced field was also discussed in connection with 
the disc's vertical structure and stability. It turns out that certain modes of instability 
can be damped whereas others are amplified by the disc's gravity [13,268]. Fully 
relativistic results employ numerical spacetimcs [191,252]. Within these solutions, the 
computation of emission-line profiles were performed; the case of moderately massive 
thin finite discs was treated by [164], and that of heavy toroids by [348].^° Another 
issue is the oscillation of the disc material excited by a small perturbation in the 
horizontal {R and </)) or vertical {z) direction [169]. In the static case, the formulae 

^ The pseudo-Newtonian approach was devised in order to introduce effects of general relativity 
into accretion models [263]. The black hole is treated as a Newtonian body whose gravitational 
field is determined by Paczyhski-Wiita potential, — M,/(r — rs). This simple expression mimics the 
gravitational field of a non-rotating BH; in particular, it reproduces the radii of marginally stable and 
marginally bound circular orbits of free test particles. Different variants of the pseudo-Newtonian 
potential were proposed in order to include rotation of the BH and to obtain better agreement with 
relativity (see [309] for references). The pseudo-Newtonian potential does not aspire to represent any 
gravitational theory, and the corresponding potential is not required to satisfy any field equations. 
This fact, however, does not diminish the practical value of the model, which applies even to regions 
with strong gravity and captures qualitative features of motion near the horizon. 

Let us mention here two examples of the effects that we observed recently in studying the properties 
of the Lemos— Letelier superposition. By generating a sequence of superpositions parameterised with 
relative disc mass, and demanding that (i) all the disc matter be interpretable as two equal counter- 
rotating streams of particles on stable time-like equatorial circular geodesies and that (ii) the inner 
disc rim be fixed at the least possible radius, it was found [310] that heavier discs of the family can 
reach considerably closer to the horizon not only in Schwarzschild r, but also in terms of physical 
measures — circumferential radius rexp(— i/j) and proper radial distance from the horizon. 
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for the frequencies of free horizontal and vertical oscillations give [311] 



11 



41/-2A 

= YZTj^ [y,R.R + '^Ry^R-(^y% + iR'^y.R), (51) 



r2i 



l-Rv,R 



^v\{l--2Ry,R) . (52) 



Without the external source, reduces to the square of the Schwarzschild orbital 
frequency, = M{M + \/B? + M^)~^ , whereas remains different due to the 
pericentre precession. It was found that perturbations in the horizontal direction 
are more important for the stability of light discs, whereas the vertical perturbations 
are more dangerous for heavy discs: for superpositions of the given type, the curve of 
marginal vertical stability is more restrictive than that of marginal horizontal stability 
for Md > 0.2296M. It even fully excludes discs with > 2M/7 [370]. Disc self- 
gravity makes oscillations of the inner parts faster and the region of the horizontal 
mode trapping somewhat smaller. 

When the black hole and the disc are allowed to rotate, manifestations of self- 
gravity could be even more pronounced, because in relativity the kinetic energy (as 
well as any other form of energy) also generates the field. Corresponding mass-energy 
currents give rise to frame dragging. Whereas the mass of the black hole is probably 
dominant in a real accretion system, the disc can bear much of the angular momentum, 
thus modifying the gravito-magnetic field of the central object significantly. This could 
be decisive for the mechanisms of rotational energy extraction. Note, however, also 
the possibility of having accreting systems where the disc/torus mass can actually be 
comparable to that of the central black hole and which can represent a transient state 
of gravitational collapse. In this context, the most popular system is a close binary 
of ultra-compact objects. Late evolutionary stages of such systems (in-spiral, merger 
and the ring-down phase) and the fate of debris have been studied, numerically and 
perturbatively, as a tentative source of gravitational waves (see [149] and references 
therein) and also as an engine behind gamma-ray bursts [86,246,282,296,298]. In 
particular, at late stages of the neutron star-black hole binary in-spiral, the neutron 
star could tidally break up [229,235,254,339,350]. 

Two-dimensional toroidal shapes similar to the one shown in figure |21 were 
originally invented mainly as a model for large and very diluted flows which probably 
exist in quasars and similar objects. These considerations have been continuously 
refined and discussed in various context (recently, e.g. [156] and references cited 
therein). As mentioned above, already in early works it was recognised that self-gravity 
imposes limits on this scheme and should be taken into account in calculations. Self- 
gravity is also relevant for a recently proposed idea [371] that non-stationary tori could 
be formed in accreting binaries. A relatively small size and steep pressure gradients 
of these configurations put them in a corner of parameter space where they may show 
high-frequency quasi-periodic oscillations [349], which makes them very attractive to 
astronomers, but formation of such structures still needs to be addressed. 

-"^^The perturbations are considered non-gravitating here, i.e. the spacetime metric is not perturbed. 
Allowing the perturbation to generate its own self-consistent field would only have a second-order 
effect on the frequencies. 
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4. Astrophysical discs in Newtonian regime 

From the point of view of gravitational theory, and also for the nature of the accretion 
flow, the most important part is the inner region near to the BH. However, the 
gravitational field of actual discs has substantial and even, possibly, dominant influence 
on their structure also at large distances (exceeding several xlO^rg; in physical units, 
rs = 2GM,/c^=2.95 x 10^{M,/Mq) cm) where the self-gravity is comparable to the 
central field. Newtonian physics is sufficient at those distances. The importance 
of the disc self-gravity is often checked by Toomre's parameter, Q = kCs/{ttT,): 
the local stability criterion with respect to axially symmetric perturbations requires 
Q > 1 [341,342] (cg is a local sound speed). Spiral structures may be produced 
already at Q ~ 2-4. Depending on the radial distribution of the thermodynamic 
quantities, global instabilities can indeed occur and play a key role in enhancing 
the accretion rate [209]. (Of course, at this point one must recall the theories of 
galactic dynamics [38].) This may facilitate angular momentum transport and provide 
a feeding mechanism for central black holes [118,141,318]. Again, we concentrate on 
discs in galactic nuclei although similar mechanisms are relevant also for discs around 
stellar-mass black holes and, to a certain degree, even for circumstellar discs in star- 
forming regions. 

4-1- Effects of self- gravity derived from simple arguments 

In addition to the classical equations of mass, momentum and energy conservation, 
the description of Newtonian gravitating fluids must include Poisson's equation, 

Ai/d(r) = A-Kp{r'), (53) 

which links the mass density distribution p in the disc and the self-gravitating potential 

{r' refers to the source and r is the field point). To understand the signatures of 
self-gravity expected in the fluid structure and motion, one can consider a purely 
rotating, steady-state system with zero viscosity and no signiflcant radial pressure 
gradient. The radial component of the Euler eqiiation (in cylindrical coordinates) is 
then: fi^i? — v, + Vr fd. In the absence of rotation, steady solutions are forbidden 
and the fluid falls towards the centre in a short time scale. Free fall can be avoided 
thanks to the centrifugal term which compensates the total gravitational attraction 
QR = —VR{vd + Vm)- If radial self-gravity is unimportant, that is if V^^^d <C Vi?z^,, 
then rotation of the fluid is determined by the BH mass only, ~ fi^ = -^Vrv,. 
This assumption underlies most disc models in astrophysics. In particular, Newtonian 
geometrically thin disc models fall within this category [280] . In the opposite situation, 
it is evident from the simple form of Euler's equation that, as the disc carries some 
mass, angular velocity is expected to increase in order to prevent accretion. This 
implies that (radially) self-gravitating discs arc expected to rotate faster than non- 
self-gravitating ones. 

From the previous estimate one can conclude that heavy (fluid) discs should 
exhibit a super-Keplcrian rotation curve (in terms of Q.{R) > 57k(-R))-^^ Keplerian 
rotation could still be maintained in a self-gravitating disc, but this would be a 
surprising result of some fine tuning among various quantities — pressure gradients, 
advection of momentum, effects of turbulent transport, etc. Neither can the possibility 

It was demonstrated [153] that the rotation curve of maser spots in NGC 1068 could be explained 
by gravitational attraction of the disc around the central BH. 
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of sub-Keplerian rotation be ruled out. Such a case could occur in the outermost disc 

regions if V/ji^d < 0; for instance due to the presence of a dense region near the edge 
or due to significant negative pressure gradients. Further, viscous stresses can switch 
the accretion. 

Let us now analyse the vertical component of the Euler equation in the same 
idealized situation as described above (again in the Newtonian limit): p~^VzP = 
— Vzf. — Vzfd. This is hydrostatic equilibrium in ^;-direction which tells us that 
central gravity acts against pressure gradients and tends to gather material around 
the equatorial plane of the disc. For a bare central point mass, the magnitude of the 
vertical component of gravity decreases with radius. As an inevitable consequence, 
the disc tends to 'flare' (its thickness being roughly proportional to the distance from 
the axis, or even steeper, and it can wildly oscillate in the vertical direction). The 
atmosphere of the disc (i.e. its less dense region away from the equatorial plane) 
gets thicker where the gravity is weaker. The geometrical semi-thickness of the disc 
then satisfies the relation ^ = ^-j^ > ^- Models show such flaring in the non-self- 
gravitating regime. Regarding the vertical direction, the disc self-gravity is significant 
if Vzt'd ^ V^i^,. For any fluid element this results in enhancement of the vertical 
gravity; thus, a vertically self-gravitating disc should be geometrically thinner and 
more dense than a non-self-gravitating one. As will be shown below, some accretion 
disc models even yield ^ < in those regions, so in these circumstances the disc 
thickness cannot exceed a certain limiting value. This is important for the disc 
irradiation (and self- irradiation) . 

Some fundamental properties of discs and tori can be understood in the apparently 
simple framework of Maclaurin spheroid theory. Early analytical investigations by 
Poincare, Dyson and, later, by Chandrasekhar [58] have been supplemented with 
recent numerical investigations (Hachisu [133,134]). These works shed light on possible 
shapes and dynamical characteristics of steady, self-gravitating polytropes in the 
inviscid limit. Both axially and non-axially symmetrical equilibrium patterns have 
been discussed. It tiirns out that equilibrium sequences can describe any conceivable 
system (depending on the polytropic index of the gas): these include multiple star 
systems, tri-axial (Jacobi) ellipsoids, some very strange shapes ('pear', 'amonit', etc.), 
but also discs, tori and rings. In the case of global solid-body rotation, discs, tori, and 
rings occupy the bottom-right part in the usual diagram of rotation frequency versus 
angular momentum. Furthermore, in a recent paper by Ansorg et al [15] a sequence 
of compressible rings is investigated with high precision. A bifurcation occurs from 
the Maclaurin sequence to the multi-ring sequence. What remains to be explored is 
the case of a full multi-ring sequence and its limit for the number of rings approaching 
infinity. An interesting connection could appear with the theory of disc fragmentation 
into individual rings [123]. The theory of Maclaurin spheroids and rings is built from 
simple arguments and is fundamental for understanding the problem of gravitating 
discs, but it still needs to be enriched with other physical ingredients. 

Disc modelling poses a difficult task even in the Newtonian regime. In addition to 
the technical embarrassment of solving the Poisson equation with sufficient accuracy, 
difficulties are connected with strong coupling between magnetohydrodynamic terms 
and radiative transfer in the medium, and with our ignorance of boundary and initial 
conditions (they cannot be sufficiently constrained by observations but are decisive 
for the evolution of flows). Furthermore, huge values of the Reynolds number 

1^ Note that even in the standard (non-self-gravitating) thin disc model, the black-body temperature 
of the disc medium does not scale proportionally to the central mass. The enormous diversity in 
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(typically > 10^^) are associated with accretion flows. In consequence, simulations 
do not currently reach the numerical resolution required to investigate simultaneously 
both the large-scale streams and the smallest swirls where kinetic energy is dissipated. 
This is another reason why (magneto-)hydrodynamical simulations of accretion flows 
should be interpreted with caution [212]. Despite the tremendous complexity of these 
problems, models and simulations have reached a fair level of credibility [140]. Thermal 
equilibrium of the disc and its possible fragmentation have been also explored [120]. 

Given the technical obstacles inherent in the problem, disc modellers often resort 
to a semi-phenomenological parameterisation of the turbulent transport by using a 
viscosity coefficient of the form i/t — aCgH, where H is the disc half-thickness and 
a < 1 is a parameter. The a-parameterisation, originally advocated in the context of 
binary accretion discs by Shakura & Sunyaev [314], remains a recipe enabling study 
of the characteristics of a mean accretion flow at the lowest level of approximation. 
However, it does not explain the process of triggering and maintaining the turbulence. 
Other viscosity laws have been proposed, for instance the /3-prescription which is 
based on analogy with turbulence observed in laboratory sheared flows and gives 
Pt oc v^R [217,295]. Also, there is the prescription of marginally stable {Q = const) 
discs with i/t cx v^H^/R, designed for weakly (locally) self-gravitating discs [28]. Very 
probably turbulent transport and its effects cannot be satisfactorily mimicked through 
simple prescriptions. Below we touch on several properties which can be derived from 
the classical a-prescription, as well as the properties oi Q = const discs. Given the 
current knowledge, one cannot decide which of these approximations is better. 

4-2. Poisson's equation in discs 

Association of the density field p{r) with a prescribed Newtonian potential is 
straightforward and always possible. Various classes of density-potential pairs can be 
formed [38,100]. Unfortunately, the reverse procedure, i.e. the inversion of Poisson's 
equation, is difhcult in general. No analytically tractable expression is available for I'd 
that would correspond to a general three-dimensional system of finite size and mass. 
Exact solutions corresponding to non-spherical systems cannot be written generally 
in a closed form if the matter distribution is two-dimensional, non-homogeneous or of 
non-trivial shape or topology. Astrophysical discs exhibit all three complications. 

As an illustration we show first, in figure (left panel), the density distribution 
in the standard disc model where self-gravity effects are not included [152]. The 
complexity of shape and density distribution is striking: discs appear to be complex 
objects for which the merit of special analytical (and especially self-similar) solutions 
is very limited or doubtful. 

An alternative approach is to consider the integral definition for the potential. 



and to use expansions of Green's function (for instance in cylindrical harmonics or 
another kind of equivalent infinite series representation). This approach is often 
employed in order to find self-similar solutions and to investigate the disc stability 
[150,196]. However, as pointed out in [70], convergence of the series is not guaranteed 

opacity, spatial and time scales, etc. does not allow for bare scaling between corresponding physical 
quantities from one system to another. As a result, Galactic BH candidates have rather different 
discs from those around super-massive BHs in galactic nuclei. Distinctions and common properties 
of these two fiducial cases of accretion flows have been reviewed by Mirabel [237] and Czerny [85]. 




(54) 
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Figure 3. The shape and density distribution of an accretion disc feeding a 
10* M0 black hole. Left: the WSG limit (cf. section lOl is indicated by ({R, z) = I 
(thick curve) . The limit is reached in the equatorial plane as close as ~ 600rs . At 
this point and beyond, the disc self-gravity is expected to play a dominant role in 
the structure and dynamics of the flow (the SSG limit and the unstable limit are 
out of range adopted in this plot). Right: The case of the SSG limit is shown. A 
standard-type steady-state model was employed [152,280] (the viscosity parameter 
is a = 0.1 and the accretion rate M = O.IMq yr~^). Lengths are plotted linearly 
and normalised to the Schwarzchild radius; the density scale is in g/cm^. 



in all points of matter distribution. Also, the expansion leads to uncertainties because 
the inferred series must be truncated. 

In order to be able to treat realistic systems, one has to seek numerical solutions. 
However, it is difficult to achieve sufhcient computational accuracy using standard 
techniques. Why is this? The first reason is that no boundary conditions are available 
in the vicinity of the disc; instead, they can be imposed only at infinity (i^d — + 
for r — > oo). Secondly, the physical boundary of the fluid is not known a priori 
and it can evolve with time. One cannot rely on the standard situation that is met 
within laboratory physics of a fluid conflned within walls, in which case the boundary 
conditions are imposed by experimental settings. The third difhculty is caused by the 
relatively large and variable aspect ratio: H/ R may be as large as unity and even 
greater at some parts of the disc, and as low as ~ 10~^ at other places (cf. figure O. 
This fact hinders discussion of discs in active galactic nuclei, where the geometrical 
thickness of the disc can vary wildly over its diameter. 

There are two practical ways of solving the problem numerically [242]. The 
first one employs straightforward determination of the potential from equation (|54|l 
everywhere in the disc. This approach can be very accurate if one is able to cope 
with integrable singularities inside the source (r — > r') and it can account for very 
distorted disc shapes and large density gradients. However, the computations involved 
are generally very time consuming. One can envisage some cell-based technique [325] 
(as done in particle simulations [143]), but this is not expected to reach high precision 
(which is necessary for stability analysis). Another possibility is to invert equation 
H53|l using a finite difference method or a finite element method [329]. This approach 
requires that boundary conditions can be obtained through equation H54|l. Finally, 
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it is worth noting that the inversion of Poisson's equation for rotating stars can be 
performed via the mapping technique, e.g. transformation of space coordinates [41,70]. 
To our knowledge, the mapping technique has never been used for discs. 

Since the efficient numerical computation of the potential is not trivial, people still 
seek approximate semi- analytical approaches that would complement purely numerical 
solutions and make them faster. This is in particular the case for the infinite slab 
(is) approximation, which was introduced by Paczyhski and collaborators in the late 
1970s in a series of papers devoted to marg inally stable discs [186, 260, 261]. In this 
approach the disc is examined in (i?, z)-plane, imposing axial symmetry and omitting 
the radial and azimuthal terms in equation H53|l . The Poisson equation then reads 

??-47rp. (55) 

The main advantage of equation H55|l is the possibility of determining the magnitude 
of the vertical gravity by direct integration along the vertical direction. In this way 
even a full solution fd(i?, z) could apparently be found. The vertical field is given by 

gti^) = -^ = -47rS(z), (56) 

where 

J:{R,z) = / p{R,z')dz' (57) 

is half of the local surface density. This result is exact if the system is invariant under 
any translation perpendicular to z-axis, and for this the disc must be (i) of infinite 
size, (ii) strictly flat (its geometrical thickness 2H is constant) and (iii) homogeneous 
in the radial direction (vertical gradients can be non-zero in the direction of z-axis). 
Obviously, such properties are not realistic, but the validity of the IS approximation can 
be extended to radially inhomogeneous, finite size discs with non-constant thickness, 
provided that the condition 

is satisfied. It is not generally possible to check inequality (|58|l for a given configuration 
because the "true" potential Vd is not known or is difflcult to compute with sufficient 
accuracy (for the technical reasons mentioned above). Also, the above-described way 
of getting i^d{R, z) cannot be used to estimate the radial component of the field. 

It is worth noting that equation H56|) holds regardless of the relative masses 
of the disc and the central object, making the IS approximation often very useful 
(for example, it can be used to compute massive non-Keplerian discs), but it is also 
important to know where this powerful approximation fails. Although the potential 
is a global function, we can determine this by considering its local properties by 
expressing z/dOci?°i |z|'^2 . Inequality lf5S|) is violated in regions where 




a,[a,-l)_^ (59) 



^"^ See also numerous follow-up papers which helped to clarify substantial uncertainties in the original 
formulation: [28,40,53,151,152,154,211,265,299,316-319,329]. In particular, the role of self-gravity 
and the disk fragmentation in fuelling the central BH has been discussed, but this topic remains 
beyond our discussion. 
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or, in terms of density gradients, when 

(Vfl, pf » (V,p + 2- V) (V.p + 2z- V) • (60) 

For a reasonable density stratification and attractive potential, we expect ai < 
and 1 < a2 < 2. This implies that the IS approximation may break down when 
(i) ai~a2, which corresponds to discs and tori of moderate and large aspect ratios, 
and/or (ii) the radial and vertical density gradients are of the same order. This applies 
if density waves or shocks with steep fronts are present, especially in boundary layers 
or close to the disc edges. The flatter is the disc, the weaker are the edge effects. 

An interesting property of vertical stratification can be derived directly from 
the IS approximation [119,260]. If the disc is able to reach hydrostatic equilibrium, 
then the total pressure gradient compensates total gravity. Differentiating equation 
of hydrostatic balance relative to z gives 



d 
di 



1 dP Mz 



This equation becomes independent of the radius in the limit of a disc that stays in 
vertical equilibrium under its own gravity. An analytical solution can be obtained 
assuming polytropic relation P — c^p^^~ (cg is a speed, n is polytropic index). 
The above-given equation is analogy of the famous hydrostatic equilibrium problem 
for gaseous spheres in equilibrium. The resulting "plane" Lane-Emden equation 
writes [123] 

Id/ 1^ dp\ ^ A-Kn ^ ^g^-j 
p dz \ dz J Cg (1 + n) 

For vertically isothermal structure with n — > oo, the equation H62|l is the Frank- 
Kamenetski equation [303] 



d^ In p Anpn 



= (63) 



dz2 c2 (1 + n) 
and the solution is [123,265,299] 

p{z) = pc cosech^ \ zc^^ ) , (64) 



where Cg is the isothermal sound speed and pc = p(0) is the mid-plane density. From 
equation (|57|l , the surface density of the self-gravitating disc is given by S = ^j2Pc/'K 
[265], and thus the disc mass depends solely on the radial distribution of central 
pressure Pc — PcCg- Discussion of solutions for various polytropic indices can be found 
in [155,260]. The issue is that, as the polytropic index n runs from infinity (in the 
isothermal case) to smaller values, the vertical density profile tends to be flatter. In 
the limit n ^ 0, this profile becomes exactly rectangular in shape. Figure0]shows the 
density profile in the direction perpendicular to the disc plane in a typical disc. The 
isothermal case is also shown for comparison. Models of the disc vertical stratification 
are close to the models of stellar interiors and show that the vertical structure is not 
isothermal and cannot be modelled by a polytrope (with n = const), at least not if 
the disc is optically thick. Discs are heated by turbulent motion and this results in a 
very complex vertical stratification of physical quantities. 
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Figure 4. The vertical distribution of density in a standard-type disc (same 
conditions as for figure |2J at distance r = 500r^ (left panel) and r = 600ts (right 
panel). Three cases are shown, all with M, = 10* Mq: (i) In the first case (thin 
solid line), hydrostatic equilibrium is maintained by the central object while the 
disc self-gravity is deliberately neglected. The mid-plane density is then equal to 
pc- These calculations include vertical heat transfer by radiation and convection, 
(ii) In the second case (thick solid line), the disc vertical self-gravity is introduced 
in the IS approximation [152]. Despite the evident accumulation of matter closer 
to the disc mid-plane, the area under the curves (which corresponds to half the 
surface density T,{H), as labelled on each curve) is almost unchanged. Details of 
stratification are sensitive to the mode of energy dissipation, (iii) Also shown is 
the solution for a self-gravitating disc in the approximation of the Lane-Emden 
equation with isothermal structure (as given by equation 1641 : dashed line). 



4-3. What criterion for self-gravity in thin discs? 

In many respects it would be useful for the disc modelling if the degree of self-gravity 
were indicated by an unmistakable criterion. As a relevant quantity, one can intuitively 
suggest the mass ratio q (— M^/AI,), where 

A/d(i?)=47r/ ^{R,H)R'dR' (65) 

is the cumulative mass of the disc and is a half of the local surface density 

according to equation l(57)l. One might postulate that self-gravity is unimportant for 
q < 10%, say. Note that, in a binary system, the q- value (defined as the mass ratio of 
the two stars) is critical in the definition of the Roche lobes. However, such a criterion 
is unreliable to explore an extended disc in its close neighbourhood, where it cannot 
be reduced to a point mass. We show below that the right criterion is based on the 
quantity ~ qR/H which contains information about the disc mass (relative to the 
central black hole) as well as about the distribution of the mass density in space (in 
terms of the aspect ratio) . 

The ratio C = g'^/g' can be evaluated provided that the vertical component of 
the disc self-gravity is known. For C = 1, the disc and the central body contribute 
comparably to the field along z-axis (this case is referred as the "weak self-gravitating 



Topical review 



29 



limit", WSG hereafter). As soon as C exceeds unity, the disc is vertically self- 
gravitating. Within the IS approximation, g'^ is given by equation (|55|l . so one obtains 
C = Att1^{R, z)r^ / {M,z). When the gas remains confined within altitudes z/R <C 1, 
we can use the mid-plane value 

Cc - -j^, (66) 

where pc is the mass density in z = 0, as defined in the absence of self-gravity. The role 
of parameter C(^) z) was discussed in [151]; Taylor expansion around z — shows that 
the mid-plane value of C can be identified with ^-parameter of Paczyhski [260] . One 
finds that the upper limit for the disc height is typically H < 0.3R; such discs stand 
on a borderline beyond which the approximation of geometrically thin discs becomes 
very crude and inadequate. We notice that vertical self-gravity occurs preferentially 
(i) when the central BH has a small mass, (ii) if the local mass density in the disc is 
large, and (iii) at rather large distances from the centre. If the surface density obeys 
a power-law in radius, for instance in the form E — Soi?", then 

Cc = g(i?) (5 + 2) (^1) {s^~2). (67) 

Relation H67|l shows that even low-mass discs can be self-gravitating if they are flat 
enough. We can put numbers in this relation: given that standard discs in active 
galaxies have a power index s ^ —0.6 and small, almost constant aspect ratio, 
H/R ~ 0.01, we find Cc = 140 x q{R) [151]. (The adopted value of s corresponds to 
the outer, gas-pressure dominated region; another value would not change qualitative 
conclusions at this point.) In other words, a disc becomes vertically self-gravitating 
when its relative mass exceeds a small threshold value. This value is model dependent; 
typically it is order of a few xlO~^. Note that the limit is much less than the crude 
estimate of 10% suggested before. 

By introducing ratio ry = V^i'd/Vfli^, one can follow the same kind of analysis 
for the radial component of gravity. A disc is radially self-gravitating if Jry] > 1, hence 
this regime can be referred to as a strongly self-gravitating (ssg) disc. Conclusions 
are less straightforward in this case because we lack sufficiently reliable expression for 
gn = — Vfd. In other words, there is no obvious analogy to the IS approximation. 
Some authors [236,346] have proposed that the total potential can be approximated 
by the monopole term, ly^ ~ —Md{R)/R, which is not fully satisfactory but gives at 
least the right order of magnitude in the SSG limit {rj ^ 1). Considering, as above, a 
power law for the radial surface density profile, rj is given by 

^ qiR) - Cc^ ^ (1 + s)q{R) (68) 

in the equatorial plane. Using equation H67|) . we find 

l + sH 

Vc = Cc 69 

(the case of s = —2 can be derived in similar manner). The SSG regime corresponds 
to 

2 + s f H\^^ , , 

Using the above values for s and H/R^ rjc = 1 corresponds to Cc ~ 150, well above 
the WSG limit. As Cc is expected to be an increasing function of radius, the SSG limit 
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involves more mass and is reached at larger radii than the weak limit. It therefore 
lies inside the vertically self-gravitating discs. This value is in good agreement with 
the results established by Lantian & Xiaoci [190] who have attempted to performed 
a direct, analytical integration of equation (|54|l . Although their investigation still 
contains some rough approximations (for instance, singularities in the Poisson kernel 
were artificially removed) , these authors derived a formula for radial self-gravity which 
accounts for finite size effects. For s = — 1, they found radial acceleration to be 



4 = 21] 



In 



{2R + H){R + Rin){Ront- R) H 



(71) 



{2R-H){R-Rin){Rout + R) 4.R\ 

which sets the SSG limit at Cc ~ 67 (assuming the Eddington accretion rate; see 
figure 01 right). Note that this formula still remains singular at the disc edge [234]. 
To circumvent the difficulty, a softening parameter is often introduced [195], however, 
correct behaviour can be obtained just by careful treatment of edges. 

Results discussed so far have not yet addressed the essential questions concerning 
the gravitational stability of self-gravitating, geometrically thin discs. It is known for 
long, and confirmed by observations, that galaxies (treated as discs of stars, gas and 
dust) are regulated by self-gravity [38,342]. Similar processes can form spirals and bar 
patterns in purely gaseous discs. It was shown in a remarkable paper by Goldreich & 
Lynden-Bell [123] that instability sets in when 

iJ>&., (72) 

where the lower limit, depends on polytropic index n (the value of ^„ = 1.1 
corresponds to a vertically isothermal disc, while ^„ = 1.75 for n 0). Note 
that this is an analogy of Toomre's original criterion, but with a slightly different 
threshold [120,128,278,341]. It is important to recaU that Goldreich & Lynden-Bell's 
criterion has been established from an analysis of growing unstable modes through a 
linear perturbation theory, within the framework of is approximation and assuming 
a non-viscous, uniformly rotating, flat disc. It is therefore expected to be somewhat 
altered when relaxing some of the above assumptions. It is easy to recognise that 
equation (|7^ also gives ^ 18-28. 

We thus conclude that instabilities are not present in a disc if its size remains 
below the weak limit. The region separating the weak, the unstable and the strong 
limit can be narrow or wide, depending on the surface density profile; or, more 
precisely, on the function Cc{R)- Again, a very nice example is provided by Saturn's 
rings which are vertically self-gravitating at their outer edge, but gravitationally stable, 
implying that the weak limit and the unstable limit do not coincide. 



4.4- The self-gravitating regime in accretion discs in active galaxies 

According to equation (|66|l . an increase of the ^-parameter is expected as the radius 
increases. This means that there is a boundary zone that separates the inner (non-self- 
gravitating) solution from the outer (self-gravitating) one. The location of this zone 
can be found without accounting for self-gravity, simply by computing C from a disc 
model [78, 128]. This approach may underestimate the magnitude of self-gravity and 
overestimate the radius of such boundary, but it provides a useful order of magnitude 
estimate. A more reliable result must include vertical self-gravity into consideration. 
This can be achieved within the IS approximation. Furthermore, in the context of 
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Q-prescription, the WSG limit is reached at Rsg=R{( — 1). For accretion rates lower 
than a certain critical value, Merit, one finds [151] 

i?sg ~ 480 a"/27M-8/27^^-2/3_2/9^-8/9^ (73) 

while 

i?,g ^ 230 a2/9m4/9M-2/3s2/3 (74) 

for accretion rates larger than Merit; here, Rsg is expressed in Schwarzschild radii, 
and the critical accretion rate is expressed in solar masses per year, Merit ~ 
0.27 QiQ^i K~^^^ fi~^^^ Mq/jt (k is opacity). Noticeable deviations are observed when 
considering temperature- and density-dependent opacities and equation of state (not 
yet specified self-consistently in equations (|73|I - H74() '). The existence of two different 
relations arises from the fact that bh discs generally possess an inner region (dominated 
by radiation pressure) and an outer (gas pressure dominated) region where the 
physical quantities (namely, the density) exhibit different variations with radius and 
the accretion rate. The case H74|l corresponds approximately to accretion luminosity 
exceeding the Eddington limit. 

Figure[Slgives a theoretical relationship between the location of wSG limit and the 
accretion rate. It was computed by different authors using the a-prescription and the 
same input parameters (relevant for accretion flows onto super-massive bhs in galactic 
nuclei). It is noticeable that no firm agreement has been achieved about the position 
of this limit, which largely fluctuates from one author to another simply because of 
different ingredients in models. Also, there arc no direct observations available to 
confirm if this wSG-transition operates in nature and if it occurs where predicted. 
However, in active galaxies there is a striking similarity between the location and 
size of the so-called broad-line region, and the value of Rsg- This correspondence 
could indicate that the clouds are formed at the distance where self-gravity becomes 
important (cf. [75] for a systematic discussion of the correlation). According to bi- 
dimensional calculations, a lower (reliable) limit appears to be Rgg ^ 530rs for a 10* 
Mq central black hole and a = 0.1. This is only 5 milliparsecs from the centre, much 
less than previously estimated (e.g. [316]). 

4-5. One dimensional models: self-similar and asymptotic solutions 

Of simple models that can be constructed analytically on the basis of a- 
parameterisation, we pay a special attention to those where self-gravity is taken into 
account in the vertical direction, while Keplerian rotation is still imposed. This 
represents even a weaker category of WSG discs, because their rotation curve is 
assumed not to be altered by the disc mass. Within this approach, the above-discussed 
infinite slab approximation has been found very convenient, although some criticism 
can be raised at this point, because these disc models typically tend to be rather 
massive. Hence they become globally self-gravitating and no more Keplerian (this 
problem occurs especially if the surface density falls faster than the square of radius). 
Paczynski's original paper [260] was devoted to the investigation of discs with constant 
ratio of the vertical self-gravity to the vertical component of the central field, i.e. 
constant value of ('-parameter (i.e. C = g^/gl, which corresponds to A in Paczynski's 
notation). The assumption of vertical self-gravity is justified by the fact that strongly 
self-gravitating discs are prone to gravitational instabilities [123, 341] which should 
make them evolve in a couple of dynamical time scales (idyn ^ fl^^). The endpoint 
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Figure 5. The location of WSG limit as a function of the accretion rate in an 
a-disc (a = 0.1). Top axis gives the logarithm of Keplerian orbital frequency (Q 
in Hz), while the bottom axis gives the corresponding value of x = logR— ^logMg 
(radius is in cm, the mass Mg is in units of 10* solar masses). The accretion rate is 
on the vertical axis. Individual curves correspond to different models as indicated. 
The WSG limit was computed without including vertical self-gravity in all models 
except Hure (1998, reference [151]; and 2000, reference [152]). The case [151] 
is one dimensional while [152] is two-dimensional. Distortions in the curves are 
given mainly by opacity which shows drastic variations with temperature. A 
reference power-law line for the standard Shakura & Syunaev (1973, [314]) non- 
self-gravitating disc is also plotted. 



of this evolution is the state of constant Toomre's Q-parameter, which is typical for 
differentially rotating discs of stars {Q can be related to C via equation (|66|l ). 

Soon it was realized that holding = const might produce some substantial 
bias [319]; a series of papers followed to elucidate the local response of the disc to 
vertical self-gravity and to study C(^) profile. Keplerian discs were investigated in 
which ^ varies as a function of radius in a self-consistent fashion (another motivation 
for this research was the chance of discovering new branches of solutions). JVIost 
explorations were carried out for steady discs, using the so-called "vertically averaged 
approximation" , where one seeks only for radial profiles of density, temperature and 
other physical quantities (without considering vertical stratification in detail). With 
additional assumptions about opacity, equation of state, pressure regime, etc., one is 
naturally lead to self-similar solutions [231,236,346]. 

As already mentioned above, vertical self-gravity makes the accreted gas more 
concentrated around the mid-plane. As a result, the disc tends to collapse in the 
vertical direction beyond some critical radius where C ^ 1 (cf. figures 01 and . This 
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/3-disc (no self-gravity) 
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17/16 
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Self-gravitating /3-disc 


-7/8 


-1/2 


-3/8 


23/8 


Q-constant disc 


-9/7 


-15/7 


6/7 





Irradiated Q-constant disc 


-5/8 


-29/16 


19/16 






Table 1. In the WSG regime, parameters of the disc are found proportional 
in the radial direction. In this table, the values of power-law index, s, are 
given for different physical quantities: mid-plane temperature Tc{R), surface 
density S(-R), geometrical thickness H{R), and self-gravity parameter ^c{R)- The 
optically thick limit was assumed and the a-, /3- and Q-constant prescriptions 
were adopted for turbulent viscosity. Self-similarity is maintained if gas pressure 
dominates in the medium with opacity coefficient and chemical composition being 
constant [28,154]. Viscosity was incorporated e.g. by [236] in a self-similar time- 
dependent model. In spite of numerous papers, it would be still premature to 
expect a general consensus regarding the viscosity model. For further details, 
see [40, 53, 119, 139, 186, 211, 299, 319, 329] . 

property is shared by all computations of geometrically thin discs, independently of 
the actual viscosity prescription (in other words, if one assumes a-parameterisation 
for viscosity, then the disc height is almost independent of a). The opposite is true 
for geometrically thick discs which tend to be thicker when self-gravity is taken into 
account [40]. If prescription for viscosity is specified, the variation with radius can 
be found in explicit form. Examples are given in table □ given the viscosity law 
and the heating mode of the disc that is on verge of self-gravitating regime, one 
finds the i?-variation of the disc thickness, temperature, surface density, etc. Note 
that ({R) increases with R in all models. This leads inevitably to the onset of self- 
gravity, subsequent reduction of the disc thickness and further consequences on many 
aspects, such as self-irradiation, dynamical stability and the emitted spectrum of the 
disc. It may also affect the disc warping [281] in regions that are shadowed from the 
central source and cannot reprocess its radiation; although weak at those distances, 
the resulting lack of heating is in favour of gravitational instability. 

Properties derived through the a-prescription deserve a comment. The R- 
dependence of the C comes out rather steep and much more pronounced than 
what is obtained for a bare standard disc (neglecting self-gravity), in which case 
^c(^) c>c i?^^/^". This means that the transition from non-self-gravitating regime to 
the fully self-gravitating one should occur in a very narrow region; the weak, unstable 
and strong limits are expected to almost coincide with each other. As a consequence 
of the steep rise of self-gravity, the disc can suffer from a runaway-type increase of 
density; one is even led to densities typical for stellar interiors and conditions suitable 
for deuterium and hydrogen burning [265] . Temperature reaches a plateau in the outer 
region of the disc plane with value [93, 128, 152,316] 

where fxmn is the mean mass per particle, m is the accretion rate normalized to 
O-IMq/yt and ao.i = ct/O.l. Some authors have argued that such uniform mid-plane 
temperatures are unphysical and they have suggested the use of a different prescription 
for viscosity [28,93]. Interestingly, this may lead back to the original idea of Paczyhski 
that viscosity is connected with the Q-parameter. 
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Figure 6. Left: time-scales in the gravitationally unstable region of an accretion 
disc [76]. One can define the mass transport time in the disc, cooling time (relevant 
for the formation of fragments), accretion time (for a star growing to a pre- 
determined mass, M* ~ IOMq), migration time (for a star sinking towards the 
centre), and free-fall time. Right: length-scales for a black hole embedded in a star 
cluster: the tidal disruption radius within which stars are disrupted, the collision 
radius within which colliding stars are disrupted, and Schwarzschild radius (i?s) 
of the central black hole. Solar-type stars (solid lines) and solar-mass red giants 
(dashed lines) have been assumed. Figure was adapted from [101, 114]. 



We conclude this section by sliowing several typical scales of the problem in 
figure These graphs are derived from basic arguments and they can help to identify 
physical mechanisms which are important from the observational point of view. Time- 
scales relevant for the gravitationally unstable regions of accretion discs are shown as 
a function of radius (in the left frame). The possibility of star formation was taken 
into account together with continuous migration and evolution in the medium (here, 
M, — 10^ Mq was assumed for the central black hole). If the stars form a dense 
(quasi-spherical) nuclear cluster in which the central bh resides, one can define the 
radius where stars are destroyed tidally or by their mutual collisions. The resulting 
value depends on M, as well as details of the stellar structure (as shown in the right 
frame). 

5. Waiting for observations and constraints 

Significant effort has been devoted to searching for black holes during the last 
few decades. Traditionally, most approaches of present-day astronomy rely on the 
electromagnetic radiation coming from cosmic sources, but there are other possibilities 
which will be surely exploited in near future (neutrinos, high-energy cosmic particles, 
and gravitational waves). The idea of detecting black holes via (electromagnetic) 
radiation of nearby luminous matter can be traced back to the celebrated lecture, 
which rev. John Michell presented in front of the Royal Society as early as in 1783.^^ 

A paper was based on that lecture and published under the title On the means of discovering the 
distance, magnitude, &c. of the fixed stars, in consequence of the diminution of the velocity of their 
light. . . in Philosophical Transactions of the Royal Society of London, 1784 74 35. 
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Figure 7. A schematic drawing of the central part of a galactic nucleus. 
Different regions of the accretion flow are shown: the innermost, radiation- 
pressure dominated region with a central BH, the intermediate region with a 
standard disc, and the outer, gravitationally unstable part. Length-scales are 
proportional to the central mass. The indicated values are valid for M. = 10* Mq 
and slightly model-dependent [75]. 



Understanding local physics of the medium is crucial if one aims to interpret complex 
and variable spectra that are observed. X rays are particularly pertinent because they 
are produced near BH, as discussed by Fabian and others [101]. Indeed, given the 
intrinsic emissivity of the medium, there are only a few factors determining the source 
spectrum on the Earth. The main factor is the black hole mass M,; the angular 
momentum of the black hole and the disc self-gravity play secondary roles. The 
geometrical optics approximation is thus adequate and one can hope to infer actual 
values of BH parameters from the motion of test matter around it. However, one has 
to be well aware of various environmental influences that hinder real observations. 

Under conditions that are thought to be typical in active galactic nuclei, the 
inner regions of accretion flows contain diluted plasma with a corona threaded by 
magnetic fields [199,279]. High-energy photons are first generated in the corona. These 
primary photons irradiate substantially colder gas of the accretion disc, where they 
are reprocessed further (figured). A fraction of photons is reflected. The temperature 
of the reflecting medium can be ~ 10'' K in galactic nuclei and about 2-3 orders of 
magnitude higher in the case of Galactic BH candidates [56,92]. The quality of photon 
reprocessing depends critically on the ionization state and other local properties of the 
medium. In this way the resulting spectrum is formed with spectral features (lines) 
superposed on top of the background continuum. X-ray spectra contain precious 
information about rapid flows from which one can infer properties of the central source 
of gravity and the state of the nearby gas. 

Manifestation of strong gravity has been searched in radiation of accretion discs 
with the main aim of finding the evidence for the central black hole [102,223,292,293]. 
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Majority of these papers deal with non-self-gravitating discs and assume that 
integrated light of the whole source is only very weakly affected by the disc gravitation 
(however, see [164,322,348]). In the case of our Galactic centre, the presence of 
standard- type (thin) accretion disc is rather improbable (e.g. Narayan [245]). Direct 
effect of the disc gravity on spectral line profiles is below current observational 
resolution for any present-day spectra of active galactic nuclei, so it can be safely 
neglected in the first approximation (we briefly discuss this effect in the following 
paragraph). It is noteworthy, however, that self-gravity of the disc has indirect 
impact on spectral features. First, as discussed in section 0] the disc geometrical 
shape is affected by self-gravity with consequences for obscuration of some parts of 
the disc and the resulting variability, especially at large values of inclination (this is 
also the case when GR effects of the central BH are maximal, so the two influences are 
mixed). Self-gravity also controls the amount of disc self-irradiation, which produces 
indirect signs in reflected light. Finally, self-gravity induces fragmentation in the 
disc medium (depending in particular on the disc cooling properties [120]) with 
fundamental consequences on the disc structure. Fragments can even lead to star 
formation or contribute to star-bursts [76,96]. The clumpy material may also result 
in much more pronounced variability observed from the source [8] . 

Let us now discuss direct effects of the combined BH-disc gravity on observed 
light rays. In the early works of Cunningham & Bardeen [83,84], a standard non- 
self-gravitating disc was assumed. These authors initiated a series of investigations 
of model spectra with the aim of revealing the imprints of GR. Since then, many 
people have advanced various modifications of the original semi-analytical approach 
[24,48,105,224,353]. Very prominent role has been played by studies of spectral 
features (lines) of iron in the X-ray band; in particular, the so called disc-line scheme 
has been widely applied, as it captures the main spectral characteristics attributed to 
strong gravity [103,369]. The line is produced by reprocessing the primary continuum 
radiation near above the inner edge of the disc, within several tens of (see [247,297] 
for recent exposition of emission mechanisms). The line intrinsic width is small, 
typically much less than the final observed profile that results from integration over 
the whole emitting surface. Standard picture of the disc thermal spectrum [257, 264] 
was reconsidered to account also for the disc gravity [192]. For detailed examination 
it would be useful to produce templates of expected line profiles or even spectral 
models that are consistent with the formation of non-axially symmetric structures in 
accretion discs as well as the presence of a black hole in the centre [138,165,348]. 
There is an exciting possibility of detecting gravitational effects on polarized light of a 
source orbiting a black hole, but this idea still awaits technology usable for its practical 
implementation [79,230]. Effects of the black hole itself on the light propagation are 
blurred in real sources in which the signal comes from different, insufficiently resolved 
regions. 

So far, only a few cases have been discovered of an accreting black hole with a line 
that attains the characteristic shape, width and variability and could be understood 
as the interplay of Doppler and gravitational shifts. The best example of such spectral 
feature originating near a bh has been identified in the nucleus of Seyfert galaxy MCG- 
6-30-15^® where the iron line was caught at the moment of high redshift [104,161]. 
Large gravitational redshift can be connected with rapid bh rotation, because the 
radius of the disc inner edge (where the line presumably originates) decreases with bh 

Cf. Morphological Catalogue of Galaxies. 
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angular momentum increasing [362]. Magnetic fields can change this straightforward 
conclusion. Also, the line is strongly variable, so its final (continuum subtracted) 
shape cannot be resolved reliably with current knowledge; but there is a good prospect 
of overpassing these deficiencies with more detailed models and using time-resolved 
spectroscopy that will be able to track rapid features [233] . Reconstruction of accretion 
disc images will then be doable using reverberation techniques [274, 327] or Doppler 
mapping [222]. This will also require good understanding of non-standard effects 
which are usually neglected, e.g. the effect of disc eccentricity [23,326]. 

As far as BH candidates in the Milky Way are concerned, the relativistic iron line 
was detected in a well-known micro-quasar GRS 1915-1-105 [225]; rapid rotation of the 
central black hole is not necessary in order to explain the data in this case. Similar 
evidence for the broad and skewed iron line has been discovered with different degree 
of credibility in about eight other sources [226] . In spite of a limited similarity of the 
relativistic iron line in these sources, we recall that recall that accretion discs have 
negligible gravity in BH candidates, and so we turn back to accretion discs galactic 
centres. 

Moving to a region upstream the accretion flow — to the distance of hundreds rg 
and farther out from the core of a galaxy, the interaction between gaseous flows and 
other astronomical bodies (stars forming a nuclear cluster) can substantially change 
the structure, accretion rate and chemical composition of the flow, so we briefly touch 
this point, too. That way, the medium feeding the BH can be enriched above the 
level of usual solar metalicity [21,64,74,77]. On the other hand, moving satellites 
feel the presence of the interstellar gas, especially in nuclei of active galaxies with 
high concentrations of stars and the relatively dense interstellar environment [31]. 
Accreting medium thus influences the motion and evolution of these bodies. Even 
though self-gravitational effects of the medium become progressively less important 
near the centre, a relatively accurate knowledge of their impact is desired in order to 
derive, reliably, BH parameters from rotational motion of the surrounding luminous 
matter. The dissipative interaction of stellar-mass satellites and the resulting orbital 
decay proceed with a time-scale that is defined by properties of both components 
of the system: in the inner disc, hydrodynamical drag dominates the orbit evolution, 
while in the self-gravitating outer region, the disc structure becomes clumpy and direct 
coUisional interaction is reduced. In even more remote parts, gravitational relaxation 
is important for the cluster evolution, so the outer cluster is dynamically distinct from 
the core. 

Syer, Clarke & Rees [332] proposed that a stellar-mass body passing the galactic 
core can be captured on a quasi-stable bound orbit around a central super-massive bh. 
In the simplest version of the model, stars have been considered as bullets crashing on 
the disc surface and passing through at hyper-sonic speed. Zurek, Siemiginowska & 
Colgate [375], Armitage, Zurek & Davies [20] and Ivanov, Igumenshchev & Novikov 
[160] examined various effects of passages (of individual stars) through diluted gaseous 
environment (slab-disc geometry was assumed; other people explored different contexts 
of this problem [64,166,255,256]). Pineault & Landry [276] discussed the statistical 
rate and distribution of impacts over the disc surface. Self-gravity of the disc medium 
was also considered [128,317,357], as well as GR effects in stellar motion [168,356], 
effects of the orbital decay of satellites by gravitational radiation and close encounters 
[137,244,287,331], or stellar coalescences [72,116,288,290]. Gradual changes of the 
physical and chemical state of satellites and of the disc itself were discussed [21]. 
Radiation torques [259] as well as star-gas hydrodynamic interaction [125, 136, 207] 




Figure 8. Left: Mean inclination I vs. eccentricity e of a satellite orbit that 
crosses a gravitating disc. Contour lines of the Kozai parameter, c = const, 
are also plotted. The orbital evolution consists of oscillations along c-contours 
and slow diffusion across them, towards the origin e = I = 0. (This example has 
eo = 0.6, /q = 70° as initial condition.) Kozai's quasi-integral is not isolating here, 
because of the orbital decay via dissipative transitions through the disc medium. 
An approximate analytic solution is also shown by a monotonic curve R, neglecting 
the influence of the disc gravity [289]. Right: Mean eccentricity e versus semi- 
major axis a of a star crossing a gravitating disc. The orbit is represented by an 
eccentricity curve with oscillations overlying the long-term decay. Also shown are 
two reference orbits having the same initial parameters (the solid curve si) or the 
same final parameters (52) as a simplified situation without the disc gravity (the 
latter two curves exhibit continuous and purely monotonic circularisation). Large 
excursions in eccentricity occur if the disc mass is non-zero, albeit it can be much 
less than the central mass (here, ~ fO~^M,). Figure adopted from [357]. 



induce different channels of angular momentum transport in the medium. The clumpy 
structure [118,189] of the outer, self-gravitating disc reduces the efficiency of direct 
star-disc collisions, while turbulence in the medium tends to strengthen dissipation 
[375]. 

Repetitive transitions bring the satellites into the plane of the disc at almost 
circular orbits [289, 332] (though, one can speculate of several effects that slow the 
process of circularisation down, so that eccentric or non-equatorial orbits may be 
maintained even at late stages [124,145,354]). Subsequent orbit evolution is reduced 
to a situation that was addressed in connection with the formation and migration of 
bodies inside proto-planetary discs [248,359,367]. The inner cluster becomes flattened 
and its structural dimension is related to the size of the disc. One can immediately 
recognise plethora of consequences connected with non-sphericity, anisotropy and 
satellite segregation in the cluster: namely, the initial mass function is modified, as 
well as the form of spectral-line profiles and stellar velocity dispersions. The onset of 
star formation leads to gradual conversion of gas into stars in the disc [76,77, 141]. 

Star-burst phenomenon has been indeed identified in some active galactic nuclei 
on ~ 10^ pc scales, which could represent outskirts of a dusty torus surrounding the 
central black hole [18,302]. Indeed, violent star formation could coincide with the outer 
parts of a massive nuclear cluster that is embedded within the dusty torus. Even if 
Md <C M,, the effect of the disc gravity on the long-term evolution of stellar orbits 
should be taken into account together with dissipation via star-disc (hydrodynamical) 
interactions. This stands in contrast with the fact that gravitation of accretion discs 
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has been neglected in most works exploring the nuclear cluster structure. The two 
influences are linked through their dependence on the disc density. A remark is worth 
making regarding the mass ratio q (which was criticised as unsuitable self-gravity 
criterion in section 14. 3|) : it is a useful measure to assess the importance of stellar 
transitions across a self-gravitating disc in galactic nuclei, because attraction of the 
disc and the effects of transitions are both increasing, roughly proportionally, with 

Circularisation of the orbits, evolution of their inclination and the resulting 
capture rate are visibly affected by the disc gravity. An important point is the range 
of time-scales: (i) dynamical period (of satellite stars revolution around the central 
mass; the shortest time-scale), (ii) period of oscillations in orbital elements (medium 
time-scale; these oscillations are due to disc gravity) and (iii) time for bringing the 
orbital plane into the disc (caused by successive interactions with the disc; long time- 
scale). The orbital decay continues also when satellites are embedded in the disc 
plane; its rate varies depending on the ability of the body to clear a gap along its 
trajectory through the medium. Different modes depend on the satellite masses, sizes 
and on the disc thickness and surface density at corresponding radius. The evolution 
of orbital parameters is an issue in extra-solar planetary research [272] , which tackles 
similar systems on much small scales — cf. the discussion of planetary formation and 
migration [19,145,208,248]. These topics have not yet been pursued to much detail 
in the context of galactic nuclear clusters embedded in the gaseous environment. 

One can understand the origin of orbital parameters oscillations in terms of 
Kozai's theory [184], taking into account the gradual adiabatic diffusion of quasi- 
integrals of motion. In order to see the reason for orbital oscillations, one borrows 
from classical mechanics' traditional approach to the disc gravitational potential Vdif) 
as a perturbation acting on the satellite motion [49,200,240]. Again, the total potential 
is supposed to be dominated by the central field of the core, while the disc contribution 
I'd can be transformed by a canonical transformation to the form which depends on 
mean orbital elements (eccentricity e, inclination /, and argument of pericentre a); the 
bar denotes mean values) and does not depend on the fast-changing variable (mean 
anomaly). See [357] for the explicit form of this new (i/j) in terms of integration (over 
the mean anomaly) along an unperturbed Keplerian ellipse; in that paper, Kuzmin's 
potential was employed as a simple model for the disc gravitational field (it would be 
quite interesting to examine a generalization in terms of Miyamoto & Nagai's potential 
[239], which can be still performed analytically). In the first order of perturbation 
theory one finds that (vd) ~ const = {vd)o', the mean semi-major axis a is also 
constant and differs from the actual osculating semi-major axis by short-period terms 
only. The mentioned averaging procedure allows us to disregard short-term effects and 
concentrate on the long-term evolution. Furthermore, owing to the axial symmetry 
of Ud{R,z) exists an additional constant of motion [184]: y/1 — cos / ~ c= const 
(figure IHl left). This way the problem of orbital evolution is reduced to the evolution of 
e and ui which are further constrained by condition (i/d) (e, u);c, a) — (fd)o- Contours 
of the perturbation potential (vd) provide a convenient representation of mean orbital 
elements in the (e, a;)-plane. Finally, with the inclusion of orbital dissipation in the 
gaseous medium, one finds that an adiabatic drift occurs and drives the orbit across 
the contours of constant {vd); subsequently, the drift provokes excursions of the orbital 
elements when a separatrix is traversed. A representative example of a trajectory is 
shown in figure |S1 (right). The interplay between Kozai's mechanism and gradual 
changes of the perturbation potential can give occasional (but large) eccentricity 
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Figure 9. Structure of a nuclear cluster in interaction with a planar, gravitating 
disc. In this model [167,330], the long-term evolution of stellar orbits is influenced 
by dissipative transitions across the disc slab and by (small) gravitational 
attraction of the disc itself. Left: Number density nt(a) of satellites within 
a shell of semi-major axis (a, a + da) (solid curve). Two populations can be 
distinguished, corresponding to the satellites of the outer cluster (dotted) and to 
the highly flattened inner cluster (dashed). Thin dotted power-law line indicates 
the equilibrium structure of the Bahcall-Wolf cluster (n oc a 

1/4). Right: The 

number density n«(a) is represented by different levels of shading here. An 
azimuthal section of the cluster is plotted. In both plots, the gradual concentration 
of satellites towards the equatorial plane and formation of a ring structure are 
evident. 

jumps; one of them is clearly visible in the figure. 

Next step would be to consider a consistent scheme for the whole cluster of 
satellites. In their classical paper, Bahcall & Wolf [22] anticipated a way to reveal 
a hypotlic'tical c;c;utral massive hole (embedded in a coUisional environment of a 
globular cluster) through stellar velocity dispersion and broadened wings of line 
profiles integrated over the bh domain of influence. In order to calculate the form 
of central cusps, the authors employed time-dependent diffusion equation. Of other 
assumptions, the most important was the spherical symmetry and isotropy of the 
cluster itself. Probability density of satellites can be expressed in terms of their semi- 
major axes: n,(a) oc a^^'^, where n*(a) da denotes the number of satellites in a shell 
with radius a and width da. Corresponding spatial density n*(r) is thus proportional 
to r~'^/^. Modifying the model assumptions (namely, that about the importance of 
collisions) leads to different but similar density profiles (power-law at large radii). For 
example, the well-known Plummer [277] sphere has a finite-density core with density 
falling off as r~^, while Jaffe [162] and Hernquist [142] models both decline oc at 
large radii. The latter case has a sound theoretical basis in the mechanics of violent 
relaxation. See also [175,270,315,368] and, for more recent discussion and references, 
see [11,333]. Further, the collisional steady-state solutions have been reconsidered as 
a model for a cusp in the center of a galaxy [99]. Conclusions of these papers can 
be compared with more advanced models that relax some of the assumptions and 
include new effects, e.g. resonances, enhanced loss-cone depletion or rotation of the 
cluster [114,173,286,291,373]. Nevertheless, self-consistent treatment of non-spherical 
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stellar systems with the effects of self-gravity and the presence of an accretion disc 
are still an open problem. The disc induced dissipative orbital decay is distinct in its 
tendency to flatten the stellar component and even form a ring of stars that could be 
resolved observationally. 

Figure El illustrates a stationary density profile of the cluster in a model where 
its members revolve around a bh and their orbits slowly decay due to dissipative 
effects. One can distinguish different populations of satellites (see the left panel): 
those which form almost spherical cluster far from the centre, outside the disc outer 
radius, and those which have been dragged into the disc plane near the centre (the 
latter form a ring). Hence it is a certain generalisation of the Bahcall-Wolf [22] cluster. 
A complementary representation is provided in the right panel of figure|2|which shows 
an azimuthal section across the cluster. Again, the inner cluster is evidently flattened 
and its dimension is tightly related to the size of the disc. It should be recalled that, 
with this approach, a relatively small volume is explored near BH, where the central 
mass has direct dynamical influence on the embedded stars and gas (in other words, 
rotational support of the stellar motion is not a negligible fraction compared to their 
random motions). The linear size of this region is a function of M, and a*, and is 
typically less than ~ 10 pc. 

Note that the log M, ~ 4 log a^, relation resembles the empirical correlation which 
has been determined in a large sample of galactic bulges [108, 121] (although the 
constant of proportionality can be different in the innermost region around the central 
BH, which is unresolved by current observations). Hence, in the inner cluster of stars, 
the obtained velocity dispersion should be determined by the joint action of gravity 
and of other agents, including dissipation in the interstellar environment. It was 
proposed [167, 330] that the effect of interaction between stars and the accretion flow 
shapes the observed M,-cr^, relation in the innermost region of some active galaxies 
(those with sufficiently dense accretion discs in the core). The model of star-disc 
interaction leads to an axially symmetric, but non-spherical central cluster, and so 
the exact value of the power-law index of the M.-a* relation depends also on the 
observer line of sight (the expected central velocity dispersion comes out somewhat 
larger for systems seen edge-on than for pole-on view) . The model therefore represents 
rather extreme case, however, it is obvious that the interaction between stars and 
the intergalactic environment should be taken into account also at larger distances 
from the core [366]. Miralda-Escude & KoUmeier [238] applied a similar idea to a 
much wider region (and hence, applicable to actual present-day observations). These 
authors propose star-disc collisions are a self-regulating process that helps to feed the 
central BH and control its growth. The infiuence of the collisions remains imprinted 
in the bulge mass-stellar velocity dispersion relation, which concerns the region in 
galactic nuclei greatly exceeding the domain of direct gravitational influence of the 
super-massive bh. 

Naively, one could expect that dissipation of satellite's orbital energy and 
momentum should always lead to a very dense nuclear cluster. However, the 
outcome may be more complicated because of non-sphericity and satellite segregation 
(caused by differing pace of orbital evolution, according to the mass function and 
the distribution of sizes of stellar bodies). For example, recently, Zier & Biermann 
have showed [374], on the basis of numerical simulations, that a flattened shape of 
the cluster may arise due purely to gravitational effect of a massive binary BH in the 
core. They found that a stellar ring is formed on a parsec scale in the orbital plane 
of the central binary, while the cluster is gradually depleted. Another approach to 
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(self-gravitating, eccentric) galactic stellar discs was applied to the double nucleus of 
M31 [300] and even the case of Milky Way was discussed recently [205]. Indeed, as 
indicated also in this paper, there is a broad variety of reasons why the compact nuclear 
cluster should acquire a non-spherical structure in the course of evolution [194,302]. 
The effects exhibited by self-gravitating accretion flows are merely one ingredient of a 
complex problem. 

6. Conclusions 

Gravitational field of black- hole discs induces great diversity of effects and gives ample 
scope for new ideas, of which we could touch only very limited section. First of all, 
as noted in the first part of this paper, even some very basic questions remain still 
open within the framework of exact general relativistic solutions for the black hole- 
disc gravitational field. Our interest in global gravitational effects connected with 
astrophysical discs has lead us to focus on systems around super-massive BHs that 
reside in nuclei of galaxies, in which case the onset of self-gravity appears inevitable. 
Recent ideas of heavy tori being formed in certain moments of neutron-star binary 
coalescence make these topics fitting also to stellar mass objects, although their "discs" 
must be very diff'erent in many respects. Because of expected instabilities, accreting 
BHs are pertinent also for forthcoming gravitational wave cxperinK^nts. 

Astrophysical black holes are supposed to be embedded in complex cosmic 
environment which consists of multiple sub-systems spanning enormous range in the 
parameter space of masses, lengths and time-scales. Gaseous discs represent just one 
of the components in this formation. Under such circumstances, analytical means 
can oSer insight into physical mechanisms which are connected with the disc, but 
specific analytical solutions lack generality and must be inevitably simplified. In 
recent decades, much effort has been devoted to clarify, which processes dominate 
the evolution and what are the applicable ranges of model parameters. It is quite 
possible that some regimes of the disc operation have been missed in previous works. 

In this respect, an attractive recent speculation suggests that intermediate mass 
black holes (with M, ~ IO^Mq) might also exist, possibly in cores of some globular 
clusters. If this new family of black holes is revealed and if they also accrete gas 
from their close environs (current observational evidence is rather scarce), it would 
provide a fascinating opportunity to study accreting black hole sources for which the 
dynamical time is of the order of ^ 1 second. In other words, the expected variability 
of such object could span over intervals that can be readily sampled over many periods 
with available technology. 

We concentrated our attention on a fairly conventional case of gas discs around 
super-massive black holes at small, sub-parsec length-scales where the motion of 
stellar-mass orbiters can serve as an accurate tool to probe the system. This 
idea requires various intervening infiuences (such as the perturbation of the central 
gravitational field caused by the disc/torus own gravity, or the orbital dissipation 
in the gaseous interstellar environment) to be understood with sufficient accuracy. 
Interesting and fundamental questions, namely, whether and where exactly the onset 
of self-gravity occurs in real accretion discs, still remain open to further research. 
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